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This international preliminary examination report has been prepared by this International Preliminary Examining 
Authority and is transmitted to the applicant according to Article 36, 



This REPORT consists of a total of 

El 



. sheets, including this cover sheet. 



This report is also accompanied by ANNEXES, i.e., sheets of the description, claims and/or drawin; 
been amended and are the basis for this report and/or sheets containing rectifications made before 
(see Rule 70.16 and Section 607 of the Administrative Instructions under the PCT). 



^^vhich have 



These annexes consist of a total of 



sheets. 



This report contains indications relating to the following items: 
Basis of the report 
Priority 

Non-establishment of opinion with regard to novelty, inventive step and industrial applicability 
Lack of unity of invention 

Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 
citations and explanations supporting such statement 

Certain documents cited 

Certain defects in the international application 

Certain observations on the international application 
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I2SI 


II 
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III 
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IV 


□ 


V 




VI 


□ 


VII 




VIII 


□ 
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Name and mailing address of the IPEA/EP 
Facsimile No. 


Authorized officer 
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INTERNATIONAL PRELIMINARY EXAMINATION REPORT 



International application No. 

PCT/FR99/02524 



I. Basis of the report 



1 . This report has been drawn on the basis of (Replacement sheets which have been furnished to the receiving Office in response to an invitation 
under Article 14 are referred to in this report as "originally filed" and are not annexed to the report since they do not contain amendments.): 

| | the international application as originally filed. 

[^[ the description, pages *"38 s as originally filed, 

pages , filed with the demand, 

pages , filed with the letter of 



pages , filed with the letter of 



the claims, Nos. , as originally filed, 

Nos. ' ^ amended under Article 19, 

Nos. , filed with the demand, 

Nos. > filed with the letter of 20 October 2000 (20.10.2000) 

Nos. , filed with the letter of 



the drawings, sheets/fig 1/12 , as originally filed, 

sheets/fig , filed with the demand, 

sheets/fig , filed with the letter of 

sheets/fig , filed with the letter of 

2. The amendments have resulted in the cancellation of: 

I I the description, pages 

[ 1 the claims, Nos. 



1 1 the drawings, sheets/fig 



3 I I This report has been established as if (some of) the amendments had not been made, since they have been considered 
— to go beyond the disclosure as filed, as indicated in the Supplemental Box (Rule 70.2(c)). 



4. Additional observations, if necessary: 
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International application No. 
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V. Reasoned statement under Article 35(2) with regard to novelty, inventive step or industrial applicability; 



citations and explanations supporting such statement 

1 . Statement 

Novelty (N) Claims 1-2 3 YES 

Claims NO 

Inventive step (IS) Claims 1-2 3 YES 

Claims NO 

Industrial applicability (IA) Claims 1 — 2 3 YES 

Claims NO 

2. Citations and explanations 



The available prior art does not disclose or suggest a 
method which takes into account the number of edges at 
each vertex to optimise the position of the resulting 
vertex, as defined in Claim 1. For this reason, Claims 1- 
23 meet the requirements of PCT Article 33(2) and (3). 
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International application No. 
PCT/FR 99/02524 



VII. Certain defects in the international application 

The following defects in the form or contents of the international application have been noted: 

Contrary to the requirements of PCT Rule 5.1(a) (ii), the 

description does not outline the relevant prior art set 

forth in document Dl (US-A-5 590 248) and does not cite 
this document. 
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Destinataire: 
VIDON, O. 

Cabinet Patrice VIDON 

Immeuble Germanium 

80 Avenue des Buttes de Coesmes 

35700 Rennes 

FRANCE 



PCT 



NOTIFICATION DE TRANSMISSION DU 
RAPPORT D'EXAMEN PRELIMINAIRE 
INTERNATIONAL 

(regie 71.1 du PCT) 



Date d'expedition 
(jour/mois/annee) 



20. 11.00 



Reference du dossier du deposant ou du mandataire 
5072.WO 


NOTIFICATION IMPORTANTE 


Demande internationale No. 
PCT/FR99/02524 


Date du d^pot international (jour/mois/annee) 
15/10/1999 


Date de priorite (jour/mois/annee) 
15/10/1998 


Deposant 

FRANCE TELECOM et al. 





1 . II est notifies au deposant que ('administration chargee de I'examen preliminaire international a etabii le rapport 
d'examen preliminaire international pour la demande internationale et le lui transmet ci-joint, accompagne, le 
cas echeant, de ces annexes. 



Une copie du present rapport et, le cas echeant, de ses annexes est transmise au Bureau international pour 
communication a tous les offices elus. 



Si tel ou tel office elu I'exige, le Bureau international etablira une traduction en langue anglaise du rapport (a 
Texclusion des annexes de celui-ci) et la transmettra aux offices interesses. 



4. RAPPEL 

Pour aborder la phase nationale aupres de chaque office elu, le deposant doit accomplir certains actes (depot 
de traduction et paiement des taxes nationales) dans le delai de 30 mois a compter de la date de priorite (ou 
plus tard pour ce qui concerne certains offices) (article 39.1) (voir aussi le rappef envpye par le Bureau 1 
international dans le formuiaire PCT/IB/301). . . 

Losrqu'une traduction de la demande internationale doit etre remise a un office elu, elle doit comporter la 
traduction de toute annexe du rapport d'examen preliminaire international. II appartient au deposant d'etablir la 
traduction en question et de la remettre directement a chaque office elu interesse. 

Pour plus de precisions en ce qui concerne les delais applicables et les exigences des offices elus, voir le 
Volume II du Guide du deposant du PCT. 



Nom et adresse postale de I'adminstration chargee de Texamen 
preliminaire international 



Office europeen des brevets 
D-80298 Munich 

Tel. +49 89 2399 - 0 Tx: 523656 epmu d 
Fax: +49 89 2399 - 4465 



Fonctionnaire autorise. 
Corcos, E 

Tel.+49 89 2399-7418 
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TRAITE DE COOPERATION EN MATIERE 

PCT 



RAPPORT D'EXAMEN PRELIMINAIRE INTERNATIONAL 

(article 36 et regie 70 du PCT) 



Reference du dossier du deposant ou du 
mandataire 

5072.WO 


voir la notification de transmission du rapport d'examen 
POUR SUITE A DONNER preliminaire international (formulaire PCT/IPEA/41 6) 


Demande internationale n° 
PCT/FR99/02524 


Date du depot international Qour/mois/annee) 
15/10/1999 


Date de priority (jour/mois/annee) 
15/10/1998 * 


Classification internationale des brevets (CIB) ou a la fois classification nationale et CIB 
G06T17/20 


Deposant 

FRANCE TELECOM et al. 



1 . Le present rapport d'examen preliminaire international, etabli par I'administaration chargee de Pexamen preliminaire 
international, est transmis au deposant conformement a Particle 36. 



2. Ce RAPPORT comprend 4 feuilles, y compris la presente feuille de couverture. 

S || est accompagne d'ANNEXES, c'est-a-dire de feuilles de la description, des revendications ou des dessins qui ont 
ete modifiees et qui servent de base au present rapport ou de feuilles contenant des rectifications faites aupres de 
I'administration chargee de I'examen preliminaire international (voir la regie 70.16 et Instruction 607 des Instructions 
admin istratives du PCT). 

Ces annexes comprennent 6 feuilles. 



3. Le present rapport contient des indications relatives aux points suivants: 

I IS Base du rapport 

II □ Priorite 

III □ Absence de formulation d'opinion quant a la nouveaute, I'activite inventive et la possibility 

d'application industrielle - .. - . - _ • - . 

IV □ Absence d'unite de I'invention 

V EI Declaration motivee selon I'article 35(2) quant a la nouveaute, I'activite inventive et la possibility 

d'application industrielle; citations et explications a I'appui de cette declaration 

VI □ Certains documents cites 

VII yE\ irregularis dans la demande internationale 

VIII □ Observations relatives a la demande internationale 



Date de presentation de la demande d'examen preliminaire 
internationale 

15/05/2000 


Date d'achevement du present rapport 

20. tl.ti 


Norn et adresse postale de I'administration chargee de 
I'examen preliminaire international: 

^ Office europeen des brevets 
D-80298 Munich 
C^ 7 Tel. +49 89 2399 - 0 Tx: 523656 epmu d 
Fax: +49 89 2399 - 4465 


Fonctionnaire autoris§ 

Viets, A (1 $ |) 

N° de telephone +49 89 2399 2577 



Formulaire PCT/I PEA/409 (feuille de couverture) (janvier 1994) 



RAPPORT D'EXAMEN 

PRELIMINAIRE INTERNATIONAL Demande internationale n° PCT/FR99/02524 



I. Base du rapport 

1 . Ce rapport a ete redige sur la base des elements ci-apres (les feuilles de remplacement qui ont ete remises a 
i'office recepteur en reponse a une invitation faite conformement a i'article 14 sont considerees dans ie present 
rapport comme "initialement deposees" et ne sont pas jointes en annexe au rapport puisqu'elles ne contiennent 
pas de modifications (regies 70. 16et 70.17).) : 

Description, pages: 

1-38 version initiale 

Revendications, 1ST: 

1 -23 regue(s) avec telecopie du 20/1 0/2000" 

Dessins, feuilles: 

1-12 version initiale 

2. En ce qui concerne la langue, tous les elements indiques ci-dessus etaient a la disposition de I'administration ou 
lui ont ete remis dans la langue dans laquelle la demande internationale a ete deposee, sauf indication contraire 
donnee sous ce point. 

Ces elements etaient a la disposition de I'administration ou lui ont ete remis dans la langue suivante: , qui est : 

□ la langue d'une traduction remise aux fins de la recherche internationale (selon la regie 23.1(b)). 

□ la langue de publication de la demande internationale (selon la regie 48.3(b)). 

□ la langue de la traduction remise aux fins de I'examen preliminaire internationale (selon la regie 55.2 ou 
55.3). 

3. En ce qui concerne les sequences de nucleotides ou d'acide amines divulguees dans la demande 
internationale (le cas echeant), I'examen preliminaire internationale a ete effectue sur la base du listage des 
sequences : 

□ contenu dans la demande internationale, sous forme ecrite. 

□ depose avec la demande internationale, sous forme dechiffrable par ordinateur. 

□ remis ulterieurement a I'administration, sous forme ecrite. 

□ remis ulterieurement a I'administration, sous forme dechiffrable par ordinateur. 

□ La declaration, selon laqueile le listage des sequences par ecrit et fourni ulterieurement ne va pas au-dela 
de la divulgation faite dans la demande telle que deposee, a ete fournie. 

□ La declaration, selon laquelle les informations enregistrees sous dechiffrable par ordinateur sont identiques a 
celles du listages des sequences Presente par ecrit, a ete fournie. 

4. Les modifications ont entraine I'annulation : 
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RAPPORT D'EXAMEN 

PRELIMINAIRE INTERNATIONAL Demande Internationale n° PCT/FR99/02524 



□ de la description, pages : 

□ des revendications, n os : 

□ des dessins, feuilles : 

5. □ Le present rapport a ete formule abstraction faite (de certaines) des modifications, qui ont ete considerees 

comme allant au-dela de I'expose de I'invention tel qu'il a ete depose, comme il est indique ci-apres (regie 
70.2(c)) : 

(Toute feuille de remplacement comportant des modifications de cette nature doit etre ind/quee au point 1 et 
annexee au present rapport) 

6. Observations complementaires, le cas echeant : 



V. Declaration motivee selon ('article 35(2) quant a la nouveaute, I'activite inventive et la possibility 
d'application industrielle; citations et explications a I'appui de cette declaration 

1 . Declaration 
Nouveaute 

Activite inventive 



Oui : 


Revendications 


1 


-23 


Non : 


Revendications 






Oui : 


Revendications 


1 


-23 


Non : 


Revendications 






Oui : 


Revendications 


1 


-23 


Non : 


Revendications 







2. Citations et explications 
voir feuille separee 
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RAPPORT D'EXAMEN Demande Internationale n° PCT/FR99/02524 

PRELIMINAIRE INTERNATIONAL - FEUILLE SEPAREE 



Concernant le point V 

Declaration motivee selon I'article 35(2) quant a la nouveaute, I'activite inventive 
et la possibility d'application industrielle; citations et explications a I'appui de 
cette declaration 

Prendre en compte le nombre d'aretes arrivant a chaque sommet, pour optimiser la 
position du sommet resultant, comme defini dans la revendication 1, n'est ni connu ni 
suggere par I'etat de la technique anterieure disponible. Pour cette raison, les 
revendications 1-23 satisfont les conditions enoncees par I'Art. 33(2) et (3) PCT. 

Concernant le point VII 

Irregularites dans la demande internationale 

Contrairement a ce qu'exige la regie 5.1 a) ii) PCT, la description n'indique pas I'etat de 
la technique anterieure pertinent expose dans le document D1:US-A-5 590 248 et ne 
cite pas ce document. 
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REVENDXCATIONS 

L Proccdd- de simplification d\)n maiilage source M form6 d'une plurali(<5 dc 
surfaces ddfinies par dcs sommcis, des faces et dcs orientations sur ces dernieres, 
Jcdit proc£de niettant cn ocuvrc une diapc dc ddcimalion par fusion dParclc, 
5 consistant h associer a une ar&c ii decimcr, d6finic par deux sommefs (11,1 2), un 
sommet unique (13), dc fafon ii obtcnir un maiilage sirnp!ifi6 M\ 
caracierisd en ce qu'il comprend une &aps dc pseudo-optimisniion npr&s Indite 
clape de decimation par fusion d'une arfites, posiltonnant ie somnicl. rcsullati! dc 
Jadite fusion cn fonction d'un crit&re tenant compte du nombre d'ar&es vives 
10 auiour dc chacun des deux somm&ts formant I'arele A fusionner, de maniirc h 
rdduirc la deviation g<5omdtrique cntre ledil maiilage source M el ledit maiilage 
siraplifie M*. 

2. Proc6d6 de simplification d'un maiilage source set on la revendicaiion l t 
caracT6ris£ en cc que ladite dtapc de pseudo-optimisation comprend une 6tape de 

15 d6nombrcmenl des arfiics vives amour des deux sommcts forraam I'aretc h 
fusionner, et une <£lape de posilionnement dudit scwunel rdsuilant, dans laquelJe 
on distingue Jes deux cas suivauts ; 

si les nombres d 7 aretes vives sont les memes auiour dcs deux 
sommcts, on place ie sommet resultant dc la fusion au milieu du 
20 segment reliant lesdits sommets (42, 44) ; 

si les nombres d* aretes vives sont difterents, on place lc sommet 
rdsultant dc la fusion sur le sommet pr£seotant Ie plus grand 
nombre d'ar^tes _ vives (43, 45). . _ 

3. ProctSde de simplification d'un. niaillage source scion J'une quelcomjue dcs 
25 rcvendi cations 1 ct 2, caracicris£ en ce qu*il comprend une dtapc dc selection 

d*une fusion d'areie h effected, paimi loutcs les fusions d'arete possibles, lenanl 
compte : 

d'au moins une information representative dc la conrburc ctefinie 
localement auiour de Tarefe considered ; 
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d'au moins une information repr£scntaiivc de la dynamique 
g£om6trique d6flnie localement. 

4. Proc6d6 dc simplification d'un maillage source selon Ja rcvcndication 3, 
caract&istS cn ce que laditc &apc de sdlection met en ocuvrc une queue dc priori^ 

5 des ar&es h fusioEincr, en fonction d'un critfere prioritaire, ladite information 
representative de la courbure, puis d'un critere sccondaire, laditc information 
representative de la dynamique g6oin£iriquc. 

5. Proc6dd de simplification d'un maillage source scion la rcvcndication 4, 
caraci6ris<S en cc que ladite dtape dc selection gcre un seuil dc courburc* scutes ]es 

10 ar&es ayant unc courbure infdrieurc audit seuil ctanl considdr&es pour 
1' application dudit crit&re sccondaire, 

ledit seuil dtant augment^ lorsqu'aucune ar6le ne pr£sente plus une courbure 
inferieure ce dernier 

6m Vroc666 dc simplification d'un maillage source scion Tunc quelconque des 
15 revendications 1 & 5 f caractdris<5 cn cc que ladite information representative dc la 
dynamique gdom&rique appariient au groupc comprenant : 
longueur de r arete considdr£c ; 

une moyenne des surfaces des faces avoisinani laditc arfiic 
consid6r<5e ; 

20 - unc moyenne des longueurs des ar6les odjaccntes aux soinmets 

« 

fonnant ladite arete consid^rfie ; 

une combinaison de longueurs d'aretes el/ou de surfaces de faces. 
7* Procdd£ de simplilJcaUon d'un niaillage source selon Tunc quelconque des 
revcndicaLions 1 & 6, canictfirisfi en ce que Ja decimation est interrompuc cn 
25 fonction d'un des crit&rcs appartenant au groupe comprenant : 

un taux de compression atteinl ; 

une complexity g6om£trique atteinte, exprimde par un nombre dc 

somniets ou de faces ; 

un seuil de courbure atteinl. 
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8. Procddd dc simplification d*un maillage source selon rune quelconque ties 
revendicalions 1 & 7, caracldrist* en ce qu'jl consiitue une diape d'initialisaiion 
d'un proc6dd d 'optimisation gdomdlrique d'un maillage. 

9. Procddc d'optiinisation g6om<Hrique d'un maillage source, caracK?/isc en 
5 ce qu'il comprend une dtape d'initialisalion mctlanl en ceuvre 1c proc6d£ de 

simplification dc Tune quelconque des revendications 1 & 7. 

10. Proc6d6 do codagc d'un maillage source (M) selon lu revendicaiion 9, 
rcpr6sentalif d'un objet cn trois dimensions, d&ivrait un maillnge simplifid (M 1 ) 
correspondant audit maillage source (M), lesdits maillages ^tani definis par un 

10 ensemble de sommets, d'aretes et/ou de faces, 

caract6ris£ en qu'il met en ceuvre uae £lape dc simplification scion Tunc 
quelconque des revendications I h 7, puis une 6tape dc minimisation d'un volume 
compris entnc lcdit maillage source (M) et ledit midllage simplifid (M*). 

11. Proc6d<5 de codage d'un maillage source scion la revendicaiion 10, 
15 caraclerisd en ce que, chacun desdits maillages 6tant defini par la position dc 

chacun de ses soznmets, Iadite eiape de minimisation assure la determination des 
positions des sommets (X u X 2 ,«.. X n ) dudii maillage simplify (M') minimisanllc 
volume V(M t M') entrc ledit maiDage source et Jedil maillage simplifies 

12. Proc<Sde de codagc d'un maillage source scion Tunc quelconque des 
20 revendications 10 et 1 1, caracL6ris6 en ce que ladite 6tapc de minimisation met cn 

ceuvre un processus itdratjf opfimisant progrcssivemcnt les positions des sommcts 
dudil maillage simplifie (MP). 

13. Proc6d£ de codagc d'un maillage source scion la revendicaiion 12, 
caractcrise en ce que ledit processus ilfimuf est interrompu lorsqu'au moins un des 

25 crilSres d 'arret sui vants est alteint : 

un nombre maximum derations ; 

une difference enlre deux vecteurs ddplacemenl succcssifs des 
positions des sommets inf6ncurc h un seuiJ pr£d<* termini (Bpsilon). 
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14. Procddd de codage d'un maillage source selon Tune quelconquc dcs 
revendicalions 11 & 13, caracteris6 esi ce que ladite dtape dc minimisation met en 
aruvre une mdthode de gradient adapt atif. 

15» Pxocdde dc codage d'un maillagc source selon la rcvendicalion 14, 
5 caracldrisd en ce que Indite mdthode de gradient adaptatif repose sur les operations 
suivantes: 

selection d'un vecteur X p de R 3n (n £ 1) dudit maillagc simplifies, c( 
calcul du gradient VE (X p ) cn X p de la fonction h minimiscr : 

E=d(M,M'(X, Xj); 

10 - determination de la position X p * de X p dudit maillage, selon la 

relation ddfinie & rildration k+1 par : 

X > " YtH ||LE(X/>i| 
k variant dc 0 h n-1 et y k dtant le pas de ladite relation. 
J6. Proc6d6 de codage d'un maillage source selon Tune quclconquc dcs 

15 revendications 10 & 15. earaclerisd en ce que, k chaque iteration, on ddlermine unc 
variation dldmentairc dudit volume conespondant fl un champ de vecteurs ° , 
el en cc que, la surface dtant paramdtrde par u et v, de faf on qu*un vecteur °- 
s*exprime sous la forme 0J * L (u, v), ladite variation dldmemaire est assimildc au 
paralldlepipcde gdn£rd par revolution de V6l6mom de surface dudv dans la 

20 direction 5] ^'(u, v). < 

17. Procddd de codage d*un maillage source selon Tunc quelconque dcs 
revendications 10 St 16* caracidrise en ce que ledil maillagc simplifid est parnni<5lr6 
h l'aide d'un module d 'elements finis. 

18* Proc&te de codage d'un maillage source selon la rcvendicalion 17, 
25 caraclcrisd cn ce que lesdits dldmenLs finis sont obtenus 2 Paide d'un interpoiatcur 
affine. 

19. Procddd de codage d'un maillage source selon I'une quclconquc des 
revendications 10 a IS, caracterisd en ce qu'il met en c&uytc un codage progrcssif 



FEUILLE MODIFI E 



20- 1 0-2000 • '• • FR 009902524 

43 



dudit maillage simplify par decimation d'un sommet et optimisation locale ou 
par ddcimaiion d'un ensemble sommets el optimisation globale. 
20. Procfiflfi de codage d'un maillage source selon I'unc queluonquc dcs 
rcvendications 10 h 19, caracterisd en oe qu'il comprend une 6tapc dc limitation 
5 de Ja degradation due k une transformation el&ncnlaire, mcltant en ocuvre une 
queue de priority sur les transformations 61dmentaires. 

2J. Proc&I6 de codage d'un maillage source selon ]a revcndicalion 20, 
carac(6risd en ce que laditc <5tapc de limitation dc la ddgradation due k unc 
transformation £l<5mentaire, c'est-il-dire une fusion d'aretc, d6finie par deux 
10 soiusneu, comprend les 6tapes de : 

calcul d'un coOt pour cbaque transformation 61£ineniaire jx)ssible ; 

realisation de la transformation dltJmentairc de moindre coflt ; 

recalcul des coflts des transformations etementaires modifies par la 

transformation cl6menlairc pr<5c£dcntc ; 
15 - ojout des nouvelles transformations 616mentaires cr66es el calcul dcs 

coQls coirespondants. 

22. Proc6d6 dc codage d'un maillage source selon la revindication 21, 
caracf£ris£ en cc que le coftt d*une transformation 6l6mentairc (Tj) s'cxprime par : 

C^CX^^^-maxd.^.FCXl)) 

20 avec : Tj transformation fusionnant deux sommets X ; et X ; dudit 

maillage simplify M* ; 
X/ le sominet dudit maillage simplifi6 M* resultant dc Indite 
transformation ; 

F(Xf)Ies faces dudit maillage simplify MVavoisinant 1c sonimel X ( ! 
25 r aprfcs ladite transformation ; 

V M ensemble des sommets dudit maillage source M appartcnant 
aux faces ayant 6l6 intersect&s lors du calcul dc r orientation 
des suifaces pendant ladite minimisation. 

23. Application du proc<5d£ de codage d'un maillage source scion Tunc 
30 qnelconque des revendications 1 h 22 & au moins un des domaines suivant : 
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la r<SaIil£ virtuclle ; 

la simulation scicntifique ; 

la moddlisation. 
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Reference du dossier du deposant ou 
du man data ire 

5072. WO 


POUR SUITE voir la notification de transmission du rapport de recherche internationale 
(formulaire PCT/IS A/220) et, le cas echeant, le point 5 ci-apres 

A DONNER 


Demande internationale n° 

PCT/FR 99/02524 


Date du depot international (jour/mois/annee) 

15/10/1999 


(Date de priorite (la plus ancienne) 
(Jour/mois/annee) 

15/10/1998 


Deposant 

FRANCE TELECOM et al . 



Le present rapport de recherche internationale, etabli par I'administration chargee de la recherche internationale, est transmis au 
deposant conformement a 1'article 18. Une copie en est transmise au Bureau international. 

Ce rapport de recherche internationale comprend £ feuilles. 

PH II est aussi accompagne d'une copie de chaque document relatif a I'etat de la technique qui y est cite. 



1. Base du rapport 

a. En ce qui concerne la langue, la recherche internationale a ete effectuee sur la base de la demande internationale dans la 
langue dans laquelle elle a ete deposee, sauf indication contraire donnee sous le meme point. 

I | la recherche internationale a ete effectuee sur la base d'une traduction de la demande internationale remise a I'administration. 

b. En ce qui concerne les sequences de nucleotides ou d'acldes amines divulguees dans la demande internationale (le cas echeant) 
la recherche internationale a ete effectuee sur la base du listage des sequences : 

| [ contenu dans la demande internationale, sous forme ecrite. 

| | deposee avec la demande internationale, sous forme dechiffrable par ordinateur. 

| | remis ulterieurement a I'administration, sous forme ecrite. 

| | ' remis ulterieurement a ('administration, sous forme dechiffrable par ordinateur. 

I I La declaration, selon laquelle le listage des sequences presente par ecrit et fourni ulterieurement ne vas pas au-dela de la 
divulgation faite dans la demande telle que deposee, a ete foumie. 

I I La declaration, selon laquelle les informations enregistrees sous forme dechiffrable par ordinateur sont identiques a celles 
du listage des sequences presente par ecrit, a ete fournie. 

2. II a ete estlme que certalnes revert dlcatlons ne pouvalent pas falre I'objet d'une recherche (voir le cadre I). 
3 - EH 11 y a absence d'unlte de I'lnventlon (voir le cadre II). 

4. En ce qui concerne le tltre, 

PH le texte est approuve tel qu'il a ete remis par le deposant. 

| | Le texte a ete etabli par I'administration et a la teneur suivante: 



5. En ce qui concerne I'abrege, 

le texte est approuve tel qu'il a ete remis par le deposant 

□ le texte (reproduit dans le cadre III) a ete etabli par I'administration conformement a la regie 38.2b). Le deposant peut 
presenter des observations a I'administration dans un delai d'un mois a compter de la date d'expedition du present rapport 



de recherche internationale. 
6. La figure des desslns a publier avec I'abrege est la Figure n 



PH suggeree par le deposant. Aucune des figures 

I | parce que le deposant n'a pas sugge>e de figure. n est a P ub,,er - 

I | parce que cette figure caracterise mieux I'invention. 
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Date du depot international 
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1 ■ [x] est notifie au deposant que le rapport de recherche internationale a ete etabli et lui est transmis ci-joint. 
Depot de modifications et d'une declaration selon 1'article 19 : 

Le deposant peut s'il le souhaite. modifier les revendications de la demande internationale (voir la regie 46): 

Quand? Le delai dans lequel les modifications doivent etre deposees est de deux mois a compter de la date de 

transmission du rapport de recherche internationale : pour plus de precisions, voir cependant les notes 
figurant sur la feuille d'accompagnement. 

Ou? Directement aupres du Bureau international de I'OMPI 

34, chemin des Colombettes 
1 21 1 Geneve 20, Suisse 
n° de telecopieur: (41-22)740.14.35 

Pour des instructions plus detaillees, voir les notes sur la feuille d'accompagnement. 

2. I i II est notifie au deposant qu ! il ne sera pas etabli de rapport de recherche internationale et la declaration a cet effet, prevue 
1 — 1 a Tarticle 17.2)a), est transmise ci-joint. 

3. I I En ce qui concerne la reserve pouvant etre formulee. conformement a la regie 40.2, a regard du paiement d'une ou 
1 — 1 de plusieurs taxes additionnelles. il est notifie au deposant que 

I I la reserve ainsi que la decision y relative ont ete transmises au Bureau international en meme temps que la requete 
1 — 1 du deposant tendant a ce que le texte de la. reserve et celui de la decision en question soient notifies aux offices 
designes. -■ . . - • 

| | la reserve n'a encore fait I'objet d'aucune decision; des qu'une decision aura ete prise, le deposant en sera avise. 

4. Mesure(s) consecutive(s) : II est rappele au deposant ce qui suit: 

Peu apres Texpiration d'un delai de 18 mois a compter de la date de priorite, la demande internationale serapubliee par le 
Bureau international. Si le deposant souhaite eviter ou differer la publication, i! doit faire parvenir au Bureau international 
une declaration de retrait de la demande internationale. ou de la revendication de priorite. conformement aux regies 
90£>/s.1 et 906/S.3. respectivement avant ('achievement de la preparation technique de la publication internationale. 

Dans un delai de 19 mois a compter de la date de priorite. le deposant doit presenter la demande d'examen preliminaire 
international s'il souhaite que I'ouverture de la phase nationale sott reportee a 30 mois a compter de la date de priorite 
(ou meme au-dela dans certains offices). 

Dans un delai de 20 mois a compter de la date de priorite, le deposant doit accomplir les demarches prescrites pour Touverture 
de la phase nationale aupres de tous les offices designes qui n ? ont pas ete elus dans la demande d'examen preliminaire 
international ou dans une election ulterieure avant rexpiration d r un delai de 19 mois a compter de la* date.de priorite ou 
qui- ne pouvaient pas etre elus parce qu'ils ne sont pas lies par le chapitre II. 
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recherche internationale 

^ Office Europeen des Brevets. P.B. 5818 Patenttaan 2 
JM NL-2280 HV Rijswijk 
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. Fax* (+31-70) .340-3016 
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NOTES RELATIVES AU FORMULAIRE PCT/ISA/220 



Les presentes notes sont destinees a donner les instructions essentielles concernant le depot de modifications selon 
('article 19. Les notes sont fondees sur les exigences du Trarte de cooperation en matiere de brevets (PCT), du reglement 
d'execution et des instructions administratives du PCT. En cas de divergence entre les presentes notes et ces exigences, ce sont 
ces demieres qui priment. Pour de plus amples renseignements, on peut aussi consulter le Guide du deposant du PCT, qui e3t une 
publication de TOM PI. 

Dans les presentes notes, les termes "article", "regie' et Instruction" renvoierrt aux dispositions du traite, de son reglement 
d'execution et des instructions administratives du PCT, respect ivement 



INSTRUCTIONS CONCERNANT LES MODIFICATIONS SELON L'ARTICLE 19 

Apres reception du rapport de recherche intemationale, le deposant a la possibilite de modifier une fois les revendications 
de la demande intern at ion ale. On notera cependarrt que, comme toutes les parties de la demande internationals (revendications, 
description et dessins) peuvent etre modifiees au cours de la procedure d'examen pretiminaire intemationaJ, il n'est generalement 
pas necessaire de deposer de modifications des revendications selon 1'article 1 9 sauf, par exemple, au cas ou le deposant so uh arte 
que ces demieres soient publiees aux fins dune protection provisoire ou a une autre raison de modifier les revendications avant 
la publication Internationale. En outre, il convient de rappeler que I'obtention d'une protection provisoire n'est possible que dan3 
certains Etats. 



Queiles parties de la demande Internationale peuvent etre modfflees? 

Selon Carta'cle 19, les revendications exclusivernent. 

Durarrt la phase intern at ton ale, les revendications peuvent aussi etre modifiees (ou modifiees a nouveau) selon 

1'article 34 aupres de ('administration chargee de I'examen pretiminaire international. La description et les dessins 

ne peuvent etre modifiees que selon I' article 34 aupres de I 'administration chargee de I'examen preliminaire international. 

Lore de f'ouverture de la phase nationale, toutes les parties de la demande internationaJe peuvent etre modifiees selon 
('article 28 ou, le cas echeant, selon Particle 41 . 



Ouand? Dans un delai de deux mois a compter de la date de transmission du rapport de recherche internationaJe ou de 1 6 mois 

a compter de la date de prior ft e, selon I'echeance la plus tardive II convient cependant de noter que les modifications 
seront re pu tee 3 avoir ete revues en temps voulu si elles parviennent au Bureau international apres I'expiration du delai 
applicable mais avant Tachevement de la preparation technique de la publication international (regie 46.1 ). 



Ou ne pas deposer les modifications? 

Les modifications ne peuvent etre deposees qu' aupres du Bureau intemationaJ; elles ne peuvent etre depose as ni 
aupres de 1'office reoepteur ni aupres de ('administration chargee de la recherche internationaJe (regie 46.2). 

Lorsqu'une demande d'examen preliminaire international a ete/est deposee, voir plus loin. 



Comment? Sort en supprimant entierement une ou plusieurs revendications, soit en ajoutant une ou plusieurs revendications 
nouvelles ou encore en modifiant le texte d'une ou de plusieurs des revendications telles que deposees. 

Une feuille de remplacement doit etre remise pour chaque feuille des revendications qui, en raison d'une ou de 
plusieurs modifications, differe de la feuille inrtialement deposee. 

Toutes les revendications figurant sur une feuille de remplacement doivent etre numerotees en chiffres arabes. Si 
une revendication est supprimee, il n'est pas obligatoire de renumeroter les autres revendications. Chaque fois que 
des revendications sont renumerotees, elles doivent I'etre de facon continue (instruction 205. b)). 

Les modifications doivent etre effectuees dans la langue dans laquelle la demande Internationale est public©. 



Quels documents dolvent/peuvent accompagner les modifications? 
Lettre (instruction 205.b)): 

Les modifications doivent etre accompagnees d'une lettre. 

La lettre ne sera pas pubfiee avec la demande international et les revendications modifiees. Elle ne doit pas etre 
confondue avec (a "declaration selon I' article 19.1)" (voir plus loin sous "Declaration selon I'article 19.1 )"). 

La lettre dolt etre redigee en anglais ou en francals, au cholx du deposant. Cependant, si la langue de la demande 
Internationale est I'anglals, la lettre doit etre redigee en anglais; si la langue de la demande InternationaJe est le 
francals, la lettre dolt etre redigee en francals. 
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NOTES RELATIVES AU FORMULAIRE PCT/ISA/220 (suite) 



La lettre doit indiquer les differences existant entre les revindications telles que deposees et les revendications telles 
que modifiees. E(le doit indiquer en particutier, pour chaque reven dication figurant dans la demande intemationaJe 
(etant entendu que des indications identiques concemant plusieurs revendications peuvent etre groupees), si 

t) la revendication n'est pas modifiee; 

ii) la revendication est supprimee; 

iii) la revendication est nouvelle; 

rv) la revendication rem place une on plusieurs revendications telles que deposees; 

v) la revendication est le result at de la division d'une revendication telle que deposee. 



Les examples suivants lltustrent la ma mere dont les modifications dot vent etre expliquees dans la lettre 
d'accompagnement: 

1 [Lorsque le nombre des revendications deposees inrtialement s'elevait a 48 et qu'a la suite d'une modification de 
certaines revendications il s'eleve a 51]: 

"Revendications 1 a 15 remplacees par les revendications modifiees portant les memes numeros; revendications 
30, 33 et 36 pas modifiees; nouvelles revendications 49 a 51 ajoutees." 

2. [Lorsque le nombre des revendications deposees inrtialement s'elevait a 1 5 et qu'a la suite d'une modification de 
toutes les revendications il s'eleve a 11]: 

Revendications 1 a 15 remplacees partes revendications modifiees 1 a 1 1 ." 

3. [Lorsque le nombre des revendications deposees inrtialement s'elevait a 1 4 et que les modifications consistent a 
suppnmer certaines revendications et a en ajouter de nouvelles]: 

"Revendications 1 a 6 et 1 4 pas modifiees; revendications 7 a 1 3 supprimees; nouvelles revendications 1 5,1 6 et 
17 ajoutees." ou 

"Revendications 7 a 13 supprimees; nouvelles revendications 15, 16 et 17 ajoutees; toutes les autres revendications 
pas modifiees/ 



4 [Lorsque plusieurs sortes de modifications sorrt faites]: 

"Revendications 1-10 pas modifiees; revendications 11 a 1 3, 18 et 19 supprimees; revendiattons 14, 15 et 1 6 
remplacees par la revendication modifiee 14; revendication 17 drvisee en revendications modifiees 15, 16 et 17; 
nouvelles revendications 20 et 21 ajoutees." 



"Declaration salon Particle 19.1)" (Regie 46.4) 

Les modifications peuvent etre accompagnees d'une declaration expliquant les modifications et precisant I' incidence 
que ces demieres peuvent avoir sur la description et sur les des sins (qui ne peuvent pas etre modifies selon 
rarticle 19.1)). 

La declaration sera pubiiee avec la demande intern at ionale et les revendications modifiees. 
Elle doit etre redigee dans la langue dans I aqua lie la demandelntematlonalo est publiee. 

Elle doit etre succincte (ne pas de passer 500 mots si elle est etabiie ou traduite en anglais). 

Effe ne doit pas etre con fondue avec la lettre expliquant les differences existant entre les revendications telles que 
deposees et les revendications telles que modifiees, et ne ta remplace pas. Elle doit figurer sur une feuille distincte et 
doit etre munie d'un titre pemnettant de I'identifier comma telle, constitue de preference des mots "Declaration selon 
rarticle 19.1)" " - - ' ■ ~ 

Elle ne doit contenir aucun commentaire denigrant relatif au rapport de recherche Internationale ou a la pertinence des 
citations que ce dernier contient. Elle ne peut se referer a des citations se rapportant a une revendication donnee et 
contenues dans le rapport de recherche Internationale qu'en relation avec une modification de cette revendication. 



Consequence du fait qu'une demande d'examen preliminaire intematlaial ait deja ete presentee 

Si, au moment du depot de modifications effectuees en vertu de I'article 19, une demande d'examen preliminaire 
international a deja ete presentee, le deposant doit de preference, lors du depot des modifications aupres du Bureau 
international, deposer egalement une copie de ces modifications aupres de ('administration charge e de I'examen 
preliminaire international (voir la regie 62.2a), premiere phrase). 



Consequence au regard de la traduction de la demande Internationale lors de I'ou vertu re de la phase natlonale 

L'attention du deposant est appelee sur le fait qu'il peut avoir a remettre aux offices designes ou elu3, lors de Touverture 
de ia phase nationale, une traduction des revendications telles que modifiees en vertu de rarticle 1 9 au lieu de la 
traduction des revendications telles que deposees ou en plus de oelle-ci. 

Pour plus de precisions sur les exigences de chaque office designe ou elu, voir le volume II du Guide du deposant 

fit ■ DOT" ' 
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5 

CLAIMS : 

1. Method for the simplification of a source mesh (M) representing a 
plurality of surfaces defined by vertices, faces and orientations of these 
faces, said method implementing a step of decimation by edge merger, 

10 consisting of the association of an edge to be decimated, defined by two 
vertices, with a single vertex so as to obtain a simplified mesh M\ 

characterized in that the method comprises a pseudo-optimizing step 
after said step of decimation by merger of a edge, positioning the vertex 
resulting from said merger so as to reduce the geometrical deviation between 

15 said source mesh M and said simplified mesh M\ 

2. Method for the simplification of a source mesh according to claim 1 , 
characterized in that said step of pseudo-optimization consists in 
enumerating the sharp edges around two vertices forming the edge to be 
merged and distinguishing the following two cases : 

20 - if the numbers of sharp edges are the same around two vertices, 

the vertex resulting from the merger is placed in the middle of the 
segment linking said vertices ; 

if the numbers of sharp edges are different, the vertex resulting 
from the merger is placed on the vertex with the greatest number 
25 of sharp edges. 

3. Method for the simplification of a source mesh according to any of 
the claims 1 and 2, characterized in that it comprises a step for the selection 
of an edge merger to be made among all the edge mergers possible, taking 
account of : 

30 - at least one piece of information representing the curvature 

defined locally around the edge considered ; 
at least one piece of information representing the geometrical 
dynamics defined locally. 

4. Method for the simplification of a source mesh according to claim 3, 
35 characterized in that said step of selection implements a queue of priorities of 
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edges to be merged as a function of a priority criterion, said information 
representing the curvature and then a secondary criterion, said information 
representing the geometrical dynamics. 

5. Method for the simplification of a source mesh according to claim 4, 
characterized in that said selection step manages a threshold of curvature, 
only the edges with a curvature below said threshold being considered for the 
application of said secondary criterion, 

said threshold being raised when there is no longer any edge having a 
curvature below this threshold. 

6. Method for the simplification of a source mesh according to any of 
the claims 1 to 5, characterized in that said information representing the 
geometrical dynamics belongs to the group comprising : 

the length of the edge considered ; 

a mean of the surfaces of the faces neighboring said edge 
considered ; 

a mean of the lengths of the edges adjacent to the vertices 
forming said edge considered ; 

a combination of the lengths of edges and/or surfaces of faces ; 

7. Method for the simplification of a source mesh according to any of 
the claims 1 to 6, characterized in that the decimation is interrupted as a 
function of one of the criteria belonging to the group comprising : 

a compression rate achieved ; 
- a geometrical complexity achieved, expressed by a number of 
vertices or faces ; 
a threshold of curvature achieved. 

8. Method for the simplification of a source mesh according to any of 
the claims 1 to 7, characterized in that it constitutes a step of initialization of a 
method of geometrical optimization of a mesh. 

9. Method for the geometrical optimization of a source mesh 
comprising a step of initialization implementing the method of simplification 
according to any of the claims 1 to 7. 

10. Method for the encoding of a source mesh (M) according to claim 
9, representing a 3D object, delivering a simplified mesh (M') corresponding 
to said source mesh (M), said meshes being defined by a set of vertices, 
edges and/or faces, 
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characterized in that it implements a step of simplification according to any of 
the claims 1 to 7, and then a step of minimization of a volume contained 
between said source mesh (M) and said simplified mesh (M'). 

1 1 . Method for the encoding of a source mesh according to claim 10 
characterized in that, each of said meshes being defined by the position of 
each of its vertices, said minimizing step provides for the determining of the 
position of the vertices (X-|, X2, X n ) of said simplified mesh (M') 
minimizing the volume V(M, M') between said source mesh and said 
simplified mesh. 

12. Method for the encoding of a source mesh according to any of the 
claims 10 and 11, characterized in that said minimizing step implements an 
iterative process progressively optimizing the positions of the vertices of said 
simplified mesh (M'). 

13. Method for the encoding of a source mesh according to claim 12, 
characterized in that said iterative process is interrupted when at least one of 
the following stopping criteria is achieved : 

a maximum number of iterations ; 

a difference between two successive shift vectors of the positions 
of the vertices that is below a predetermined threshold (e). 

14. Method for the encoding of a source mesh according to any of the 
claims 1 1 to 1 3, characterized in that said step of minimization implements 
an adaptive gradient method. 

15. Method for the encoding of a source mesh according to claim 14, 
characterized in that said adaptive gradient method relies on the following 
operations: 

the selection of a vector X p of R 3n (n > 1 ) of said simplified mesh 
and the computation of the gradient VE (X p ) in X p of the function 
to be minimized E = d(M,M f (X-|,... X n )): 
- the determining of the position of X p * of X p of said mesh according 
to the relationship defined in the iteration k+1 by : 

v *♦«_* " v VE(X p k) 
Xp =X * - YkX ||VE(X p k )|| 

k varying from 0 to n-1 (with n < N) and Yk being the step of said 
relationship. 
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16. Method for the encoding of a source mesh according to any of the 
claims 10 to 15 characterized in that, at each iteration, an elementary 
variation of said volume corresponding to a vector field SM is determined and 
in that, since the surface is parametrized by u and v so that a vector SM \s 
expressed in the form SM (u, v), said elementary variation is likened to the 
parallelepiped generated by the evolution of the surface area element dudv in 
the direction SM ( u , v). 

17. Method for the encoding of a source mesh according to any of the 
claims 10 to 16, characterized in that said simplified mesh is parametrized by 
means of a model of finite elements. 

18. Method for the encoding of a source mesh according to claim 17, 
characterized in that said finite elements are advantageously obtained by 
means of a refined interpolator. 

19. Method for the encoding of a source mesh according to any of the 
claims 10 to 18, characterized in that said method implements a progressive 
encoding of said simplified mesh by decimation and local optimization. 

20. Method for the encoding of a source mesh according to any of the 
claims 10 to 19, characterized in that it comprises a step of limitation of the 
deterioration due to an elementary conversion implementing a priority queue 
on the elementary conversions. 

21 . Method for the encoding of a source mesh according to claim 20, 
characterized in that said step of limitation of the deterioration due to an 
elementary conversion, namely an edge merger, defied by two vertices, 
comprises the steps of : 

computing a cost for each possible elementary conversion ; 
carrying out the lowest cost elementary conversion ; 
recomputing the costs of the elementary conversions modified by 
the previous elementary conversion ; 

adding the new elementary conversions created and computing 
the corresponding costs. 

22. Method for the encoding of a source mesh according to claim 21, 
wherein the cost of an elementary conversion (Tj) is expressed by : 

CCT^X^)) = max d 2 (V M ,F(Xf)> 
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With : 

Tj conversion merging two vertices Xj and Xj of said simplified mesh 
M' ; 

Xjf the vertex of said simplified mesh M' resulting from said 
conversion ; 

F(Xjf) the faces of said simplified mesh M' neighboring the vertex Xjf 
after said conversion ; 

Vm set of the vertices of said source mesh M belonging to the faces 
having been intersected during the computation of the orientation of 
the surfaces during said minimization. 

23. Application of the method for the encoding of a source mesh 
according to any of the claims 1 to 22 to at least to the following fields : 
virtual reality ; 
scientific simulation ; 
modelling. 
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(57) Abstract 



The invention concerns a method for simplifying a mesh source M us- 
ing a decimation step by edge fusion, and comprising a pseudo-optimisation 
step following said decimation step by ridge fusion, positioning the vertex 
resulting from said fusion so as to reduce the geometric deviation between 
said source mesh M and said simplified mesh M\ Advantageously, said step 
consists in counting the sharp edges around two vertices forming the edge to 
be merged, and in distinguishing between the two following situations: if the 
number of sharp edges are the same around two vertices, the vertex resulting 
from the fusion is set in the middle of the segment joining said vertices (42, 
44); if the number of edges are different, the vertex resulting from the fusion 
is set on the vertex presenting the highest number of sharp edges (43, 45). 

(57) Abrege' 

L'invention conceme un proc&te de simplification d'un maillage 
source M, mettant en oeuvre une 6tape de decimation par fusion d'arete, 
et comprenant une etape de pseudo-optimisation apres ladite 6tape de 
decimation par fusion d'une arete, positionnant le sommet resultant de ladite 
fusion de maniere a reduire la deviation geom6trique entre ledit maillage 
source M et ledit maillage simplifid M* avantageusement, cette &ape consiste 
a denombrer les aretes vives autour des deux sommets formant Tarete a 
fusionner, et a distinguer les deux cas suivants: si les nombres d'aietes vives 
sont les memes autour des deux sommets, on place le sommet resultant de la 
fusion au milieu du segment reliant lesdits sommets (42, 44); si les nombres 
d'anites vives sont differents, on place le sommet resultant de la fusion sur 
le sommet presentant le plus grand nombre d'aretes vives (43, 45). 
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Procede de codage d'un maillage source, avec optimisation de la 
position d'un sommet resultant d'une fusion d 'arete, et applications 
correspondantes. 

/ . Domaine de ^invention et applications 

LI Domaine de V invention 

Le domaine de 1' invention est celui du codage de structures de donn^es 
geom&riques, ou maillages, notamment de grande taille. Plus precisement, 
I'invention concerne la representation et le codage d'objets ou de scenes en trois 
dimensions. Plus precisement encore, I'invention concerne une technique 
d'approximation d'un maillage source en trois dimensions pouvant etre 
utilisee seul, ou en combinaison avec d'autres techniques connues. Dans ce 
dernier cas, le procede de I'invention peut constituer une etape avantageuse 
d'initialisation. 

Un maillage est classiquement defini par un ensemble de sommets et de 
faces orientees definissant une topologie. De tels maillages sont par exemple 
utilises en graphisme sur ordinateur, pour modeliser des objets en trois dimensions 
avec une complexity geom&rique limitee. 

L' approximation d'un maillage M consiste k trouver un maillage M' 
dont la complexite geometrique est moindre que celle du maillage M, et qui 
approche au mieux la geometrie de M. 
12 Exemples d' applications 

L'invention trouve des applications dans tous les domaines ou il est 
souhaitable de reduire le nombre d'informations necessaires pour representer et/ou 
manipuler efficacement un objet en trois dimensions ou un ensemble d'objets, par 
exemple pour 1' analyser, le stocker et/ou le transmettre et/ou en assurer le rendu. 

A titre indicatif, I'invention peut notamment s'appliquer au domaine de : 
- la halite virtuelle (visites ou boutiques virtuelies, loisirs, 

telemanipulation, etc.). Dans ce type duplication, l'approximation de 

maillages permet de reduire le coflt de rendu de scenes complexes, en 

particulier en definissant la notion d'Schelonnabilite sur les maillages 
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(fonction du point de vue, des capacites graphiques, du taux de 
rafraTchissement souhaite, ..'.). Dans le cas de la realite virtuelle 
distribute ou partagee, cela permet egalement d'adapter la complexite 
d'une scene aux capacites de rendu et de stockage des differents 
5 terminaux, ainsi qu'aux debits des reseaux ; 

- la simulation scientifique (elements finis, CAO, etc.). La reduction de 
la complexite geome"trique des modeles permet une acceleration des 
temps de calcul, une prise de decisions plus rapide, notamment lors de 
la conception en CAO, et 1'elimination des informations redondantes 

10 dans une base de donnees 3D ; 

- la moderation (scanner 3D (reconstruction de surfaces a partir de 
points non organises), scanners volumiques, reconstruction de surfaces 
a partir de photos stereoscopiques ou de sequences video, de modeles 
numeriques de terrains (imagerie satellite ou radar), etc.). Un modele 

15 numerique de terrain permet ainsi l'obtention d'un maillage 

representant la topologie d'une region. Un tel maillage est obtenu par 
l'echantillonnage regulier d'une image stockant Information d'altitude 
en chaque point. II en r&ulte une quantite de donnees importante, 
comprenant des informations inutiles pour la simulation scientifique, ou 
trop couteuse pour le rendu (dans le cas de simulateurs). 
L'approximation de maillages rtduit la quantite de donnees, tout en 
garantissant une bonne fid^lite geomttrique aux donnees initiates et la 
conservation de la topologie. 
2.Artanterieur 

25 2J. Les families d'al^orithmes 

Plusieurs techniques d'approximation de maillages sont deja connues. 
Les plus repandues peuvent etre classees en trois grandes families d'algorithmes 
selon qu'elles fonctionnent par : 

decimation ; 

30 - re-echantillonnage sous-critique ; 
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subdivision adaptative. 
2.7.7 decimatinn 

La decimation consiste a retirer de maniere iterative des sommets et/ou 
des faces d'un maillage. Cette operation est appelee operation elementaire de 
simplification. Les m<Sthodes mettant en ceuvre ce principe de decimation peuvent 
element optimiser les positions.des sommets apres ou pendant la simplification, 
cette derniere etant choisie de maniere a preserver au mieux la topologie du 
maillage. 

2.7.2 re'-echantilln nnafte snus-critig up 

Le re-echantillonnage consiste a echantillonner un modele original, soit 
en prenant des points aleatoirement sur sa surface et en retriangulant ensuite soit 
en definissant une grille tridimensionnelle et en agglomerant les sommets dans 
chaque boite elementaire de la grille. Le modele ainsi genere est simplifie, et doit 
approximer au mieux les donnees initiales. Cette technique est rapide, mais ne 
conserve pas la topologie ni les caracteristiques visuellement importantes des 
maillages. 

2 .7 J subdivision adaptativp. 

La subdivision adaptative commence avec un modele comportant une 
geometrie tres simple, que l'on subdivise recursivement ensuite, en ajoutant a 
chaque iteration un detail dans les regions ou l'erreur d'approximation est 
maximale. 

2.1.4 combinations d'alfinritft m ^ 
De facon a permettre une approximation d'un maillage avec une quality 
de reconstruction satisfaisante, il est necessaire de combiner une decimation et 
25 une optimisation des positions des sommets conserves. En d'autres termes, 
1'objectif de base d'un procede de codage d'un maillage source etant de 
maximiser la qualite de 1' approximation pour une complexity geometrique 
donnfe, celui-ci doit notamment presenter les proprietes suivantes : 

decimation ; 

30 - preservation de la topologie ; 
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optimisation des positions, suivant un critere d'en-eur 
predefini. 

On connait ainsi une premiere methode appetee « remaillage », 
repondant a ces criteres. Elle est notamment presentee dans le document « re- 
tiling polygonal surfaces » par Greg TURK (SIGGRAPH 92 Conference 
proceedings, p. 55-64, 92). Elle fonctionne par 6chantillonnage, decimation et 
optimisation des positions. Elle permet de parameter le nombre de niveaux de 
resolution souhaites et le nombre final de sommets pour chacun d'eux. 

L'algorithme debute par une distribution aleatoire du nombre de 
sommets paramdtrds sur la surface du modele, puis optimise leur position en les 
liant par des forces de repulsion (fonctions de la courbure locale). II triangule 
ensuite les polygones ainsi formes, en incluant les sommets du maillage original. 
Enfin, on decime le modele un sommet apres Tautre (si cette operation conserve 
la topologie locale) puis on retriangule le trou r&ultant. 

On obtient ainsi un modele raffine dans les regions de forte courbure et 
decime sur les regions planes. 

Un inconvenient de cette^ m&hode est qu'elle pr^sente une partie 
aleatoire. En consequence, deux decimations successives ne donnent pas le meme 
r&ultat. Par ailleurs, cette technique a tendance a lisser les discontinues dans la 
plupart des situations. 

De plus, la fidelite aux donnees initiales n'est pas exprimee directement 
au cours de T optimisation des positions. 

Une autre technique, appelfe "codage progressif de maillages", a 6t6 
developpee par Hugues HOPPE, dans le document "Progressive meshes" 
(SIGGRAPH 96 Conference Proceedings, pages 99-108, 1996). Elle repose sur la 
decimation et Toptimisation des points. Cette technique est decrite plus en detail 
en annexe 1 , en relation avec les figures 24 & 27, de fa9on k ne pas surcharger la 
pr&ente discussion. II est clair cependant que cette annexe 1, ainsi que Tannexe 2, 
font partie intdgrante de la description. 
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Cette technique de codage progressif de maillages presente differents 
inconvenients. Notarament, elle ne conserve pas naturellement les discontinues 
ou singularites geom^triques, ce qui suppose une gestion de cas particuliers, tels 
que les coins, les aretes vives finissantes et les aretes vives regulieres, pour eviter 
5 des ruptures de topologie. Cela entraine des calculs complexes, difficiles a mettre 
en ceuvre, et suppose des param&rages peu ais<5s (qu'est-ce qu'une arete vive ? 
Comment ponderer le terme de ressort pour eviter de lisser les surfaces ? ...). 

Par ailleurs, elle n'est pas efficace pour caracteriser les differences 
perceptuelles entre deux maillages. 
10 3. Obiectifs de V invention 

L'invention concerne en particulier la technique de decimation mise en 
ceuvre, par exemple, par ces differentes techniques de codage de maillage, ou par 
d'autres. 

L'invention a notamment pour objectif de pallier les differents 
15 inconvenients des techniques connues. 

Ainsi, un objectif de l'invention est de fournir un precede" de simplification 
de maillage par decimation (fusion d'arStes), qui soit plus efficace, en termes.de 
qualite perceptuelle, que les techniques connues. 

En d'autres termes, un objectif de l'invention est de fournir un tel precede 
presentant, jusqu'a un niveau eleve" de decimation, une conservation des 
singularites des maillages, et une preservation de la topologie. 

L'invention a egalement pour objectif de fournir un tel procede\ qui soit 
simple a mettre en ceuvre, en termes de calculs a effectuer, et qui presente une 
bonne vitesse d 'execution. 

Selon un premier aspect de l'invention, un autre objectif est de fournir un 
tel procede, qui puisse Stre utilise seul, afin de fournir une mdthode rapide de 
simplification de maillages. 

Selon un deuxieme aspect de l'invention, un objectif est de fournir un tel 
procede" de decimation, qui puisse etre utilise, de fa ? on a l'am&iorer, dans un 
procedd d'optimisation geom&rique de maillage. 
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Selon encore un autre aspect de 1'invention, un objectif de 1'invention 
est de foumir un precede de codage d'un maillage source en trois dimensions, 
prtsentant un meilleur rapport qualite de rapproximation/complexite 
geometrique que les techniques connues. 

Notamment, 1'invention a pour objectif de fournir un tel precede, 
respectant plus fidelement les caractenstiques importantes a conserver 
(singularites ou discontinuites). 

Un autre objectif de 1'invention est de fournir un tel precede, ne 
necessjtant pas une recherche et une gestion des cas particuliers (tels que les 
aretes et les coins), ni la mise en ceuvre de parametres speafiques (par exemple 
criteres de reconnaissance d'une arete, seuils, ...) qu'il est necessaire de definir 
pour chaque maillage. 

L'invention a encore pour objectif de fournir un tel precede, qui permet 
une reconstruction successive du maillage, celui-ci etant reconnaissable, sous 
une representation grossiere, tres rapidement. 

L'invention a egalement pour objectif de fournir un tel precede 
permettant de repondre aux fonctionnalites fixees dans le projet de 
standardisation ISO-MPEG4, a savoir : 

- la compression de maillages, qui doit reduire fortement la quantite de 
donnees decrivant la geometrie, avec des pertes visuelles minimales ; 

- la representation progressive d'un maillage, permettant de deTmir un 
maillage sous la forme d'un maillage de base a faible resolution, et 
d'une sequence de raffinements ; 

- l'adaptation a une gamme de debits importante, la qualite etant 
acceptable, mais pour de tres faibles debits. 

4. caracteristiquex pri ncipals de V inventing 

£1 — tiude d'opti misation de la position des sommets 

Ces objectifs, ainsi que d'autres qui apparaitront plus clairement par la 
suite, sont atteints selon 1'invention a l'aide d'un procede de simplification d'un 
maillage source (M) representatif d'un objet en trois dimensions, delivrant un 
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maillage simplifie (M') correspondant audit maillage source (M), lesdits 
maillages etant definis par un ensemble de sommets, d'aretes et/ou de faces, ledit 
procede mettant en oeuvre une etape de decimation par fusion d'arete, consistant a 
associer a une arete h decimer, dgfinie par deux sommets, un sommet unique, de 
fapon a obtenir un maillage simplifie M'. 

Ce procede comprend une etape d'optimisation, ou de pseudo-optimisation, 
apres ladite etape de decimation par fusion d'une arete, assurant le positionnement 
du sommet resultant de ladite fusion de maniere a reduire la deviation 
geometrique entre ledit maillage source M et ledit maillage simplifie M\ 

On respecte ainsi les particularites, notamment les aretes vives du 
maillage. 

Ladite etape de pseudo-optimisation peut avantageusement consister a 
denombrer les aretes vives autour des deux sommets formant l'arete a fusionner, et 
a distinguer les deux cas suivants : 

si les nombres d'aretes vives sont les memes autour des deux 
sommets, on place le sommet resultant de la fusion au milieu du 
segment reliant lesdits sommets ; 

si les nombres d'aretes vives sont differents, on place le sommet 
resultant de la fusion sur le sommet presentant le plus grand 
nombre d'aretes vives. 
42 caracteristiques avantaeeuses de V invention 

L'invention peut notamment etre mise en oeuvre a l'aide d'un algorithme 
organise en deux parties : l'optimisation de positions et la decimation. 

La decimation permet d'obtenir un maillage geometriquement simplifie, k 
partir du maillage original, tout en conservant la topologie et une bonne 
ressemblance avec ce maillage original. 

L'optimisation des positions des sommets permet Tobtention de la 
meilleure approximation au sens de la metrique basee volume de l'invention. 

En combinant decimation et optimisation, on peut generer soit un maillage 
sous forme progressive (en utilisant alternativement la decimation d'un sommet et 
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une optimisation locale), soit des niveaux de resolution distincts (en utilisant 
successivement la decimation d'un ensemble de sommets et une optimisation 
globale). 

De facon avantageuse, le precede de l'invention comprend une etape de 
selection d'une fusion d'arete a effectuer, parmi toutes les fusions d'arete 
possibles, tenant compte : 

d'au moins une information representative de la courbure definie 
localement autour de 1'arSte consideree ; 
- d'au moins une information representative de la dynamique 

geometrique definie localement. 
La prise en compte de ces deux criteres permet, comme on le verra par la 
suite, d'optimiser le choix des fusions a effectuer, en supprimant en priority les 
elements les moins significatifs perceptuellement. 

De facbn preTerentielle, ladite etape de selection met en oeuvre une queue 
de priorite des aretes a fusionner, en fonction d'un critere prioritaire, ladite 
information representative de la courbure, puis d'un critere secondare, ladite 
information representative de la dynamique geometrique - 

Cette hierarchie de criteres permet d'atteindre une bonne efficacite. 
Avantageusement, ladite etape de selection gere un seuil de courbure, 
seules les aretes ayant une courbure inferieure audit seuil etant considers pour 
Implication dudit critere secondaire, ledit seuil etant augmente lorsqu'aucune 
arete ne presente plus une courbure inferieure a ce dernier. 

Selon differents modes de realisation particuliers, ladite information 
representative de la dynamique geometrique peut appartenir au groupe 
comprenant : 

longueur de 1'arSte consideree ; 

une moyenne des surfaces des faces avoisinant ladite arete 
consideree ; 

une moyenne des longueurs des aretes adjacentes aux sommets 
formant ladite arete consideree ; 
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une combinaison de longueurs d'aretes et/ou de surfaces de faces ; 
toute autre grandeur caracteristique reliee a la densite locale. 
Comme on le verra par la suite, la prise en compte de la longueur de l'arete 
est une technique simple, et fournissant de tres bons r&ultats. 
5 La decimation peut notamment etre interrompue en fonction d'un des 

cri teres appartenant au groupe comprenant : 

un taux de compression atteint ; 

une complexity geometrique atteinte, exprimee par un nombre de 
sommets ou de faces ; 
10 - un seuil de courbure atteint. 

Selon un premier mode de mise en oeuvre de I'invention, le precede de 
simplification d'un maillage source constitue une etape d'initialisation d'un 
precede d 'optimisation geometrique d'un maillage. 

L'invention concerne egalement un tel precede d'optimisation geometrique 
15 d'un maillage source, comprenant une etape d'initialisation mettant en oeuvre le 
proc£d£ de simplification decrit ci-dessus. 

Selon un second mode de mise en oeuvre de l'invention, le precede de 
simplification d'un maillage source peut etre utilise seul. 

Ainsi, le precede de codage de maillage peut mettre en met en oeuvre une 
20 etape de minimisation d'un volume compris entre ledit maillage source (M) et 
ledit maillage simplifie (M'). 

Ainsi, I'invention concerne une technique d'approximation de maillages 
utilisant une metrique basee volume, et non comme cela est classiquement le cas , 
une metrique considerant la distance entre un sommet et une surface (une telle 
25 distance n'etant d'ailleurs pas unique). 

Cette approche nouvelle prend en compte naturellement les 
caracteristiques importantes a conserver (singularites) sans qu'il soit necessaire de 
recourir a la detection de cas particuliers. Aucun parametrage particulier n'est 
requis, avant le traitement d'un maillage : on minimise simplement le volume 
30 entre les deux maillages. 
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De facon avantageuse, chacun desdits maillages etant defini par la position 
de chacun de ses sommets et des aretes les reliant, ladite etape de minimisation 
assure la determination des positions des sommets (X„ X 2 ,... X„) dudit maillage 
simplify (M') minimisant le volume V(M, M') entre ledit maillage source et ledit 
5 maillage simplify. 

Prtferentiellement, ladite etape de minimisation met en ceuvre un 
processus iteratif optimisant progressivement les positions des sommets dudit 
maillage simplifie" (M'). 

Ledit processus iteratif peut notamment gtre interrompu lorsqu'au moins 
10 un des criteres d'arret suivants est atteint : 

un nombre (N) maximum d'iterations ; 
- une difference entre deux vecteurs deplacement successifs des 
positions des sommets inferieure a un seuil predetermine (e). 
Selon un mode de realisation avantageux, ladite etape de minimisation met 
en oeuvre une methode de gradient adaptatif. Cette mdthode, ainsi qu'on le verra 
par la suite, permet de mettre en ceuvre Invention avec une complexity de calcul 
reduite. En effet, il n'est alors pas necessaire d'effectuer un calcul explicite du 
volume. 

Dans ce cas, ladite methode de gradient adaptatif peut avantageusement 
20 reposer sur les operations suivantes: 

- selection d'un vecteur X p de R 3 " (n a 1) dudit maillage simplify et 

calcul du gradient VE (X„) en X,, de la fonction a minimiser 
E=d(MM'); 

- determination de la position X/ de X„ dudit maillage, selon la relation 
25 definie a 1 'iteration k+1 par : 

Y k+i v k VE(X k ) 



15 



|vb(x/: 

k variant de 0 a n-1 (avec n < N) et Yk etant le pas de ladite relation ; 

N reprdsente le nombre d'iterations utiles pour atteindre la convergence (cette 

derniere etant atteinte lorsque la difference entre deux determinations 
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successives est inferieure a un seuil ( £ ) ou lorsqu'un nombre maximum 
d' iterations est atteint): 

On peut noter que X p appartient a R 3 si on effectue une optimisation locale 
ou a R 3n si on effectue une optimisation globale. En d'autres termes, on optimise 
les positions d'un sommet a la fois dans le cas d'une optimisation locale et de tous 
les sommets simultanement (sous la forme d'un unique vecteur de R 3n ) dans le cas 
d'une optimisation globale. 

De facon prtferentielle, ledit pas Vk est variable, et il varie en fonction de 
I'oscillation de deux deplacements successifs des sommets et/ou en fonction des 
variations de l'energie. II peut notamment 6tre maintenu entre deux bornes et 

Ymax" 

Selon un mode de realisation avantageux de Invention, on determine, a 
chaque iteration, une variation elementaire dudit volume correspondant a un 
champ de vecteurs ^ . La surface etant parametree par u et v, de facon qu'un 
15 vecteur «M s ' e xprime sous la forme ^(u, v), ladite variation elementaire est 
assimilee au parallelepipede genere par revolution de l'element de surface dudv 
dans la direction 8M " (u, v). 

Cette variation elementaire de volume entre M et M' peut 
avantageusement s'ecrire : 

5d(M,lvr ,5lvr ) = -/ J t](u,v)n(u,v)5M (u,v)do(u,v) 

n(a,v) = r,(M,Nf ,n(u,v)) - 1 , a fi(u>v) ^ 

orientee vers le volume interieur, et -1 sinon. 
Selon un mode de realisation preferentiel de l'invention, ledit mailiage 
simplify est param6tre a l'aide d'un modele d'elements finis. Lesdits elements 
25 finis sont avantageusement obtenus a l'aide d'un interpolateur affine. 
Dans ce cas, ledit mailiage simplifie peut s'ecrire • 

N 

M(u,v)=£X i X i (u,v) 



20 



i-l 



avec X ; (u,v) une fonction de forme adaptee a un modele d'elements 
finis triangulaires, 
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et Xj sommet dudit maillage, defini dans R 3 , 

et on definit une fonction de forme sur le maillage par le biais de coordonnees 
barycentriques. 

En utilisant Texpression prec6dente, on obtient : 

N 

6d(M,M ,ax, 6X D ) — // ri(u,v)n(u,v) £ SX^vWu.v) . 



u,v 

N 



i-l 



~ - E J I rj(u,v)n(u,v)8X i X i (u,v)da(u,v) 

im{ u,v 

De fa?on avantageuse, la position d'un sommet X { a la k* mc iteration 

s'ecrit : 

dbd 



ax : 



avec la derivee partielle de la distance pour un sommet Xj : 
dbd 

-rrr - " ! ! ri(u,v)n(u,v)X i (u,v)da(u,v) 

u.veSuppflO 

Avantageusement, le proe^de de T invention met en oeuvre un codage 
progressif dudit maillage simplifie, par decimation et optimisation locale. Ainsi, 
ledit maillage simplifie est, de fa?on avantageuse, represents par un maillage de 
base et une sequence de raffinements dudit maillage de base. 

De fa9on avantageuse, le calcul du gradient a la k ifemc iteration comprend 

les Stapes suivantes : 

- discretisation de l'expression de la derivee partielle de la distance pour 

chaque sommet X i? sous la forme : 
dbd 11 

"av""~E E I T l( 1 xJx) R ( i xJx)^i( i TJT)da(i T ,jJ 

reS i -1 j -1 

T X 

avec : - S Tensemble des triangles voisins du sommet X i; 

N T le nombre de points echantillonnes dans la direction de u ; 

le nombre de points echantillonnes dans la direction de v ; 
calcul de 1'orientation des surfaces, par identification de la 
plus proche intersection avec le maillage source M. 
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Selon un aspect preferentiel de Tinvention, l'orientation relative des 
surfaces audit maillage source et dudit maillage simplifie est calculee selon 
1'equation : 



T|(u,v) = -<n M .,n M >.<n M ., XjX M > 

avec : Xj le point echantillonne sur M' ; 

X M le point d'intersection de la droite passant par X } et de 

direction n avec le maillage source M ; 

la normale au maillage source M au point X M ; 

la normale au maillage source M' au point Xj ; 
X X 

10 ( j m represente un vecteur). 

Selon un autre aspect avantageux de l'invention, le precede comprend une 
etape de limitation de la degradation due a une transformation elementaire* 
mettant en oeuvre une queue de priorite sur les transformations elementaires. 

Preferentiellement, ladite etape de limitation de la degradation due a une 
15 transformation elementaire comprend les etapes de : 

calcul d'un cout pour chaque transformation elementaire possible ; 
realisation de la transformation elementaire de moindre cout ; 
recalcul des couts des transformations elementaires modifiees par 
la transformation elementaire precedente ; 
20 - ajout des nouvelles transformations elementaires creees et calcul 

des couts correspondants. 
Le cout d'une transformation elementaire (T { ) peut par exemple s'exprimer 

par : 

CCT^X^Xp) = max d 2 (V M ,F(Xf)) 
25 avec : T { transformation fusionnant deux sommets X t et Xj dudit 

maillage simplifie M' ; 
Xi f le sommet dudit maillage simplifie NT resultant de ladite 
transformation ; 

FCX/) les faces dudit maillage simplifie M' avoisinant le 
30 sommet X/ apres ladite transformation ; 
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V M ensemble des sommets dudit maillage source M 
appartenant aux faces ayant ete intersectees lors du calcul 
de Torientation des surfaces pendant ladite minimisation. 
Le procede de Tinvention peut s'appliquer k de tres nombreux domaines 
5 techniques, et en particulier a au moins un des domaines appartenant au groupe 
comprenant : 

la reality virtuelle ; 
la simulation scientifique ; 
la moderation. 
10 5. liste des figures 

D'autres caracteristiques et avantages de Tinvention apparaitront plus 
clairement a la lecture de la description suivante d'un mode de realisation 
preferentiel de Tinvention, donnee a titre de simples exemples illustratifs et non 
limitatifs, et des dessins annexes parmi lesquels : 
15 - la figure 1 illustre le principe d'une fusion d'arete ; 

la figure 2 illustre le principe de la queue de priorite combinant la 
courbure et la dynamique geometrique definies localement, selon 
Tinvention ; 

la figure 3 illustre le denombrage des aretes vives autour d'un 
sommet du maillage ; 

la figure 4 presente le principe de la pseudo optimisation entre fusion 
d'aretes, selon Tinvention. 

la figure 5 illustre le principe, consistant a minimiser le volume entre 
deux surfaces ; 

les figures 6 a 11 sont discutees en annexe 2, decrivant les aspects 
mathematiques de Tinvention, et illustrent respectivement : 

- figure 6 : variation elementaire du maillage ; 

- figure 7 : variation elementaire du maillage dans la direction de 
la normale ; 

- figure 8 : systeme de coordonnees sur un triangle ; 
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- figure 9 : fonction de forme ; 

- figures 10 et II : exemples d'e volutions de surface sur deux 
courbes en deux dimensions ; 

la figure 12 illustre la mise a jour des adjacences de face apres fusion 
d'une arete adjacente h. deux faces ; 

la figure 13 illustre la mise a jour des adjacences de face apres fusion 
d'une arete adjacente a une face ; 

la figure 14 illustre Techantillonnage sur un triangle pour 
Tintegration numerique ; 

la figure 15 explique le principe de T orientation relative des 

surfaces, de fa$on a eviter la formation de plis sur la surface lorsque 

la position initiale est trop eloignee de T optimum ; 

la figure 16 illustre le cout d'une transformation elementaire ; 

la figure 17 est un synoptique de Tapproximation de maillages selon 

Tinvention ; 

la figure 18 est un synoptique de T optimisation d'un sommet Xi ; 

les figures 19 a 22 illustrent le comportement du procede de 
Tinvention, sur un parallelepipede rectangle, apres decimation sur un 
coin (figure 19), une arete vive reguliere (figures 20 et 21) et un plan 
(figure 22) ; 

la figure 23 illustre les differentes (12) iterations effectuees sur un 
exemple simple, selon Tinvention ; 

les figures 24 a 27 sont discutees en annexe 1, presentant la 
technique connue du codage progressif de maillages, et illustrent 
respectivement : 

- figure 24 (egalement commentee par la suite) : transformation 
elementaire (fusion d'arete) ; 

- figure 25 : solutions initiates ; 

- figure 26 : mise a jour des transformations dans le voisinage 
affecte par la derniere modification ; 
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- figure 27 : ajout des transformations nouvellement creees. 
6, principe general d' une mise en ceuvre de Vinvention 
6.0 Introduction 

Selon l'aspect principal de Tinvention, on prevoit une £tape d'optimisation 
5 de la position des sommets (§ 6 r 2), qui Sanscrit dans le cadre d'une fusion d'arete, 
choisie par exemple selon la technique decrite au § 6.1. Cela peut ensuite etre mis 
en ceuvre dans le cadre des procedes des § 6.3 et 6.4. 
6.L Simplificati on eeometrig ue 

Ainsi qu'indique precedemment, 1'invention concerne notamment une 
10 nouvelle technique de simplification d'un maillage 3D, reposant sur la mise en 
ceuvre d'une queue de priorite combinant la courbure locale et la dynamique 
geometrique locale. Cette technique presente notamment Tavantage de conserver 
jusqu'a un niveau eleve de decimation, les singularites sur les maillages. EUe 
possede de plus une vitesse d'execution interessante. 
15 Selon un aspect de Tinvention, on construit une queue de priorite g&ant 

Toperateur topoiogique de fusion d'aretes. Cette queue de priorite combine les 
criteres de courbure locale et de dynamique geometrique locale, afin d 'exploiter 
au mieux le degre de liberte donne par Tordre des transformations a realiser sur le 
maillage. 

20 La simplification d'un maillage M consiste a construire un maillage M' de 

complexity geometrique reduite, qui conserve une faible deviation geometrique 
avecM. 

L'algorithme de simplification geometrique doit permettre de specifier une 
resolution geometrique au sommet pret. Pour cela, on choisi un operateur 
25 topoiogique Slementaire de simplification presentant de bonnes proprietes : 
conservation de la topologie dans une certaine mesure de decimation, et absence 
de creation de trous sur les surfaces, et conservation des orientations. 

Get operateur topoiogique elementaire est la fusion d'aretes, tel que defini 
par exemple par HOPPE (document deja cite) qui est illustre en figure 1 . 
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La fusion d'aretes 10 consiste h. fusionner les deux sommets adjacents 1 1 
et 12 en un sommet 13, k supprimer les deux faces 14 et 15 et a positionner le 
sommet 13 resultant de la fusion. On notera que cette transformation est reversible 
(possibility d'insertion 16 d'un sommet). 

Chaque transformation 61ementaire decime le maillage approximant M\ 
La qualite de Fapproximation se degrade done au cours de la decimation, ou reste 
au mieux invariants Afin de limiter les degradations apport^es au maillage, on 
souhaite bien sur effectuer tout d'abord les transformations affectant le moins 
possible le modele. 

Pour cela, on definit une queue de priorite contenant toutes les 
transformations realisables sur le maillage (soit approximativement le nombre 
d'aretes). Durant la decimation, la transformation de moindre cout de la queue de 
priorite est realisee, puis supprimee de la queue. Le cout dans le voisinage modifie 
par r operation precedente est alors recalcule, et les nouvelles transformations 
potentielles sur le maillage sont inserees dans la queue de priorite, apres en avoir 
calcule le cout. 

Les singularites geometriques, qui sont des parties tres informatives, 
doivent etre conservees le plus longtemps possible pendant la decimation. En 
particulier, les regions de forte courbure apparaissent comme tres informatives. En 
consequence, le premier critere de tri sur les transformations elementaires est done 
lie a la courbure locale autour de Tarete a fusionner. 

On appelle courbure C(X) autour d'un sommet X; l'angle maximal entre 
les normales a deux faces adjacentes autour du sommet X 4 . On appelle ensuite 
courbure autour d'une arete (reperee par deux sommets X ; et Xj) la moyenne de 
ces criteres evaluee en chaque sommet. 

Un second critere base sur la longueur de l'arete a fusionner, permettant de 
reduire la densite geometrique du maillage et d'obtenir un bon rapport d'aspect 
sur les triangles resultants, est mis en ceuvre. 

On peut egalement utiliser une formule usuelle representative de la 
compact. Cependant, la longueur de l'arete a fusionner pr^sente les avantages : 
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de fournir un maillage de densite uniforme dans les regions de 
courbure voisine ; 

de conserver une bonne compacite des triangles, puisque ce critere 
tend a creer des triangles equilateraux ; 
5 - de presenter un faible coutdecalcul. 

En d'autres termes, la queue de priorite selon I'invention repose 
simultanement sur la prise en compte des deux aspects suivants : 

un petit triangle n'est interessant que dans une region hautement 
informative (soit une region de forte courbure) ; 
10 - i] est souhaitable de reduire la density d'un maillage afin d'en 

reduire sa complexite. 
Les deux criteres pris en compte selon I'invention sont combines de 
maniere a obtenir le comportement illustre par la figure 2. Cette figure est une 
echelle de la courbure, graduee de 0 a n, en radians. 
15 Sur * es regions de faible courbure 21 , inferieures au premier seuil 22, la 

densite est reduite, et la compacite obtenue est raisonnable, puisque les aretes de 
longueur minimale sont fusionnees. 

Le seuil definissant une faible courbure est ensuite augmente (23), lorsqu ! il 
n'y a plus de transformation possible sur le segment de faible courbure 2L 
20 Ainsi > sur les niveaux de decimation les plus faibles, la contrainte de 

courbure est automatiquement relachee, afin d'atteindre la complexite 
geometrique fixee. 

La queue de priorite possede done deux niveaux de contrainte, organises 
suivant une hierarchie : la courbure en est le critere prioritaire, et la densit6 en est 
25 le critere secondaire. 

62 Pseudo-optimisation 

Apres une operation de fusion d'aretes, on s'autorise, selon Taspect 
essentiel de Tinvention, a placer le sommet resultant de la fusion sur la position de 
Toptimum probable. 
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Pour cela, on introduit la notion d'arete vive : une arete est vive lorsque 
Tangle forme par les normales aux deux faces adjacentes est superieur h un seuil 
fixe parametrable. On denombre ensuite le nombre d'aretes vives autour des 
sommets de Tarete & fusionner, ainsi que cela est illustre en figure 3. Le sommet 
5 31 n'est assocte a aucune arete vive, le sommet 32 a deux aretes vives et le 
sommet 33 a trois aretes vives. 

On deduit de ce denombrement deux cas pour 1'initialisation, ainsi que 
cela est illustre en figure 4 : 

- si les nombres d'aretes vives autour de X a et autour de XfeSont 
identiques, on effectue Initialisation au milieu du segment 
formant Tarete a fusionner. Sur les zones planes de faible courbure, 
cela permet de conserver une bonne compacite sur les triangles 
voisins. Sur une arete vive reguliere (nombres d'aretes vives egaux 
a 2), cela permet de positionner le sommet proche de Tarete vive 
qui sera conserve au cours de Toptimisation ; 
si les nombres d'aretes vives autour de X a et autour X b sont 
differents, Initialisation est effectude sur le sommet presentant le 
plus grand nombre d'aretes vives. Dans les cas les plus courants, 
Toptimum est atteint a partir de cette position initiale. 
Dans Texemple de la figure 4, dans lequel le maillage source correspond a 
un parallelepipede 41 , on constate que cette heuristique place le sommet : 

sur le coin du parallelepipede lorsque Tarete a fusionner forme un 
coin 43 ; 

et sur Tarete vive reguliere du parallelepipede lorsque Tarete 
commence sur la region plane et se termine sur Tarete vive 45. 
Les situations ou le nombre d'aretes vives est le meme autour des deux 
sommets sont illustres en 42 et 44. 
63 Applications 
Comme indique precedemment, la technique de simplification de 
Tinvention peut etre mise en oeuvre seule, pour offrir une technique 
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d'approximation de maillage, ou comme une etape d'initialisation d'une 
procedure plus complete d'optimisation geometrique, telle que par exemple celle 
decrite ci-apres. 

63 J structure 

5 ^invention peut etre mis en ceuvre dans une methode d'approximation de 

maillage mettant en ceuvre une m6trique basee volume. Un des avantages de cette 
metrique est qu'elle prend naturellement en compte les caracteristiques 
importantes a conserver (les singularites), sans devoir recourir a la detection de 
cas particuliers. 

10 D'une fa?on generate, 1' approximation de maillages genere un modele 

polygonal geometriquement plus simple que le modele original. Le nombre de 
polygones necessaire a la representation est ainsi reduit, tout en conservant la 
meilleure approximation, au sens perceptuel, du modele original. 

Le probleme correspondant peut etre formule par une approche 
15 variationnelle permettant d'approximer un maillage M par un maillage M' (ou % 
comprenant un nombre reduit de triangles. Elle autorise le deplacement des 
sommets en optimisant leurs positions au sens d'un ou plusieurs criteres definis 
par le biais d'une fonctionnelle d'energie a minimiser. 

L' approximation d'un maillage genere un modele geometriquement plus 
20 simple que le modele original. Le nombre de polygones necessaires a la 
representation est ainsi reduit tout en conservant la meilleure approximation (au 
sens perceptuel) du modele original. La qualite de T approximation definit le 
niveau de ressemblance perceptuelle pour une complexite geometrique fixee. 
L'invention utilise une metrique basee volume entre "le maillage approximant" et 
25 le "maillage approxime". 

Le principe de l'algorithme est scinde en deux parties : la decimation et 
1'optimisation des positions. 

L'operation elementaire de decimation choisie est la fusion d'arete, decrite 
precedemment. L'algorithme d'optimisation choisi est un algorithme iteratif de 
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gradient a pas adaptatif afin de minimiser it^rativement le volume compris entre le 
modele et le maillage original. 

En combinant decimation et optimisation, on peut generer soit un maillage 
sous forme progressive, soit des niveaux de resolution distincts. 
5 632 decimation 

Afin de simplifier la geometrie d'un maillage, il convient de decimer 
iterativement le maillage tout en conservant une bonne topologie. On choisit done 
une transformation Elementaire : la fusion d'arete, qui fusionne deux sommets 
adjacents en un sommet, supprime deux faces et place le nouveau sommet 
10 resultant de la fusion, suivant le principe decrit precedemment. 

Cette decimation peut etre rfialisee alternativement ou successivement avec 
une optimisation (telle que decrite precedemment) des positions d'un ou de 
plusieurs sommets suivant F application visee. On peut construire plusieurs 
niveaux de resolution d'un maillage en utilisant cette transformation elementaire, 
15 et la granularite minimale est obtenue par une decimation puisque le maillage est 
alors defini au sommet pres. On procede alors en alternant decimation locale et 
optimisation locale du sommet resultant de la fusion. 

Inversement, on peut definir des niveaux de resolutions de maniere 
generate en decimant au prealable un ensemble de sommets puis en optimisant 
20 globalement les positions des sommets du maillage. 

La decimation successive permet de construire des maillages progressifs. 
Dans ce cas, il convient de choisir judicieusement l'ordre de realisation des 
transformations. Cet ordre est defini par le biais d'une queue de priorite sur les 
transformations elementaires de maniere a effectuer en premier les 
25 transformations affectant le moins possible le modele original. 

Cela amene a definir la notion de cout pour chaque operation elementaire 
et la queue de priorite sera triee par ce dernier. Ainsi initialement, toutes les 
transformations possibles sont listees, le cout de chacune d'elles est calcule (sans 
effectuer de modifications sur le maillage), on realise la transformation de 
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moindre cout, et on recalcule les cofits des transformations elementaires modifiees 
par la transformation precddente. 

Bien sGr, d'autres techniques que la decimation peuvent etre utilises pour 
obtenir le maillage approche, ou approxime, tel que des techniques de sous- 
echantillonnage, de subdivisions adaptatives, etc... 

optimisation basee sur le volume 

Le probleme de roptimisation consiste done a minimiser une fonctionnelle 
d'erreur definie par : 

E(M,M') = d(M,M') 
oil d(M, M') caracterise une distance entre les deux maillages M et M\ et garantit 
la fidelite aux donnees initiales. L'algorithme d' optimisation est choisi notamment 
pour ses performances en terme de vitesse et/ou de pertinence de la solution 
finale. 

On explique ci-dessous le contexte theorique de la metrique basee volume 
utilisee pour caracteriser 1'erreur d'approximation entre deux maillages, ainsi que 
la formulation variationnelle permettant 1'optimisation de la qualite de 
l'approximation, selon 1'invention. 

6.33.1 metriq up 

Un exemple de mesure de 1'erreur, ou distance, entre les maillages M et 
M' est illustree en figure 5. La mesure de 1'erreur est donnee par le volume V, 
compris entre les deux surfaces formees par les maillages Met M'. 
Ce volume est deTmi par l'integrale de Lebesgue suivante : 

V(M,M , )=/.I v (M',M)(q)d 3 q 

25 avec v(9) Ja fonction indicatrice de rint&ieur du volume et q un point 
deR 3 . 

6.33.2 fonctionnelle. 
Un maillage 6tant defini partiellement par les positions de ses sommets, on 
peut formuler le probleme de la minimisation du volume par le biais d'une 
30 fonctionnelle d'energie a minimiser, cette fonctionnelle 6tant dgfmie par : 



20 



WO 00/22577 



23 



PCT/FR99/02524 



E = E mr = d(M, M') = d(M,M\X lv ..,X D )) 
L'algorithme d 'optimisation mis en ceuvre consiste alors a determiner les 
positions optimales des sommets X, , X 2 , ...X n qui minimisent la distance 
d(M\M). 

6333 variation elementaire de volume 
Involution elementaire du maillage M' est definie par le champ de 
vecteurs SM? . Une expression de la variation elementaire du volume est decrite en 
annexe 2 (§ 2). 

63 3 A mise en ceuvre du gradient adaptatif 
La technique d'optimisation utilisee pour resoudre le probleme de la 
minimisation peut avantageusement reposer sur une approche iterative, tel que le 
gradient adaptatif. L'annexe 2 (§ 1) decrit en detail la mise en ceuvre de cette 
technique du gradient adaptatif, ainsi que le calcul d'une variation elementaire de 
volume permettant la determination du gradient a chaque iteration, et le 
parametrage du maillage par un modele d'elements finis triangulares. 

Un modele d'elements finis triangulares permettant de calculer le gradient 
en chaque sommet est decrit en annexe 2 (§ 3). 

L'algorithme du gradient pour un sommet X; s'ecrit sous la forme : 
v k + i V k dSd 

L'expression de la derivee partielle de la distance pour un sommet X ; est : 

~7Z~- J J v)n(u,v)X i (u;v)do(u 1 v) 

u,v eSupp(\) 

Le support de la fonction de forme est defini par les triangles voisins du 
sommet Xi. L'integrale est calculee par le biais d'une discretisation sur les 
triangles (decrite par la suite) et une heuristique a ete adoptee pour determiner 
Torientation relative des surfaces formees par M et M' (egalement decrite par la 
suite). 

A chaque iteration du gradient, on applique un deplacement sur le sommet 
Xi qui tend a minimiser le volume entre les surfaces. Le deplacement est la 
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somme des contributions liees aux calculs d'orientation en plusieurs points de 
chaque triangle. 

Le pas Yi est adaptatif. Au depart il est initialise" a Yo (une fraction de la 
longueur de l'arete fusionnSe), puis il est multiplie par k (avec k < 1), en presence 
d'une oscillation (c'est-a-dire lorsque deux vecteurs de deplacement successifs 

dbd k " 35d k 

<_ T Y , y > < 0 

sont de directions opposees : OA ' dA i ). 

Ce pas ^ est borne" (borne inferieure Y min ), afin de conserver la convergence 
de ralgorithme. On a Yi ]5, 2^6], avec 0 < 6 < ^.et HHE^IIslAo (ou HE 
represente le Hessien de E). 

On considere la convergence atteinte lorsque le deplacement du sommet X; 
d'une iteration a l'autre est inferieure a e (eetant une fraction de la dynamique 
locale des aretes du maillage). Le nombre d'iterations de ralgorithme peut 
egalement etre borne. 

6 3.4. description d'urj mode dp. realisation dAtailU 
Le precede" de 1'invention permettant de generer des approximations de 
maillages a differents niveaux de resolution, un codage progressif de maillages 
par decimation et optimisation locale peut etre mis en ceuvre. 

Dans cette optique, on definit une transformation elementaire de 
decimation adaptee au codage progressif, une methode d 'optimisation locale et 
une queue de priorite sur les transformations afin de maximiser le ratio quantity 
d'informations/qualite visuelle a chaque iteration de ralgorithme. On decrit par la 
suite la mise en ceuvre de cet algorithme, et notamment le calcul du gradient deTmi 
sur un sommet du maillage; 

63.4.1 tran sformation elementaire 
Le codage progressif permet de definir la notion d'echelonnabilite sur un 
maillage. L'echelonnabilite maximale est obtenue lorsque la resolution du 
maillage peut etre choisi au sommet pres. Un maillage peut ainsi etre code sous la 
forme d'un maillage de base (a une resolution visuellement acceptable) et d'une 
sequence de raffinements. 
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Afin d'obtenir une telle structure de donnees, il convient de decimer 
iterativement le maillage, d'optimiser localement les positions des sommets de 
maniere a approximer localement au mieux le maillage initial, et d'enregistrer 
simultanement la sequence de decimation afin de coder le maillage dans le sens 
desraffinements. 

On choisit done une transformation elementaire reversible, par exemple la 
fusion d'aretes d^finie par HOPPE (document dej& cite), qui est illustree en figure 
1. La fusion d'arete 10 consiste h fusionner les deux sommets adjacents 11 et 12 
en un sommet 13, k supprimer les deux faces 14 et 15 et a optimiser la position du 
sommet 13 resultant de la fusion. 

Cette transformation est reversible (insertion 16 d'un sommet). 

Apres transformation, les adjacences doivent etre mises a jour. Pour cela, 
on considere une structure de type faces/sommets et une adjacence definissant 
toute les relations entre les elements. Les faces stockent des liens vers leurs faces 
voisines, et chaque sommet liste les faces qui le contiennent et les sommets 
voisins. Les figures 12 et 13 illustrent la mise a jour des adjacences de faces apres 
fusion d'une arete adjacente a deux faces (figure 12) ou a une face (figure 13). 

La decimation est ainsi naturellement prise en compte sur toutes les zones 
du maillage, sans qu'il soit necessaire de detecter les cas particuliers, telles que les 
aretes adjacentes a une seule face. 

6342 initialisation 

L'optimisation agit sur la position du sommet resultant de la fusion de 
r arete formee par les sommets X, et X b . Afin de debuter l'algorithme 
d'optimisation avec une bonne condition initiale, on choisit de positionner le 
sommet le plus proche de Toptimum probable. Pour cela, on introduit la notion 
d'arete vive : une arete est vive lorsque Tangle forme par les normales aux deux 
faces adjacentes est superieur a un seuil fixe parametrable. On denombre ensuite 
le nombre d'aretes vives autour des sommets de l'arete a fusionner, selon 
l'approche deja discutee. 

6343 optimisation 
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On rappelle que l'expression de la fonctionnelle consideree : 

E^E^d^M') 
Cette fonctionnelle est minimisee en optimisant la position du sommet X; 
par un algorithme de gradient : 

xr=xr- Yi — 

5 dX i 

Les determinations du pas et des criteres d'arret sont discutees au §- 

6.33.4. 

63 4.4 discretisation 

Le calcul du gradient necessite done une discretisation de F integrate : 
35d 

— — J J T](u,v)n(u,v)X i (u,v)da(u,v) 

10 ° A i u^eSupp^ 

Sur le support de la fonction de forme definie sur Xj. 

On est done amene a echantillonner la surface des triangles, avant 

d'evaluer en chaque point X(x,y,z) la valeur de la fonction de forme ^(x, y, z) et 

r orientation des surfaces ti(u,v). 

15 Le gradient s'exprime done sous la forme discrete suivante : 

38d T x 

"a3T=-Z I E ^(ixJx^axj^ia^jjdaa.jj 

U/V i xeSi-lj-l 

t J x 

avec : - S l'ensemble des triangles voisins du sommet X { ; 

- N, le nombre de points echantillonnes dans la direction de u ; 

- le nombre de points echantillonnes dans la direction de v. 

20 Pour un triangle donn£, Fechantillonnage est effectue dans le plan passant 

par ses sommets, le nombre de points etant proportionnel a la surface du triangle, 
ainsi que cela est illustre en figure 14. 

Le repere d'echantillonnage est forme par l'arete 161 de plus grande 
longueur et la hauteur 162 du triangle courant. Soit n le nombre minimum de 
25 points par triangles fixes a Favance, S min l'aire minimale des triangles voisins de 
Xj et h le pas d'echantillonnage. On deduit h par la formule suivante : 

nxh 2 =S n 
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La fonction de forme est evaluee par un rapport de surfaces tel que decrit 
en annexe 2. 

63.4 5 orientation des surfaces 
L' orientation relative de deux surfaces definie en trois dimensions est 
delicate & definir avec exactitude. On utilise avantageusement une heuristique 
pour calculer I'orientation des deux surfaces M et M' (c'est-a-dire le terme r \( u ^ v ) 
de Tequation du gradient presente ci-dessus) : 

ti(u,v) - - < n*,fi M > . < n w , X } X M > 

avec : Xj le point echantillonne sur M' ; 

X M le point d' intersection de la droite passant par Xj et de direction 
njvr avec le maillage source M ; 

nM la normale au maillage source M au point X M ; 

n ™ la normale au maillage source M' au point Xj. 
ou <,> represente le produit seal aire. 

La figure 15 illustre ce calcul de Torientation en deux dimensions. En trois 
dimensions, il suffit de remplacer les segments par des triangles. Cette heuristique 
limite fortement la formation de plis sur la surface, puisque le terme <n M , ' n M > 
traduit Fopposition des normales lorsqu'elles sont de sens opposes. 

63.4.6 priorites sur les transformations 
Chaque transformation elementaire decime le maillage approximant M'. 
La qualite de Tapproximation se degrade done au cours de la decimation, ou reste 
au mieux invariante lorsque le volume reste inchange, et ce malgre roptimisation 
des points. 

Par exemple, si Ton considere une sphere, on sait qu'elle necessite en 
theorie une infinite de triangles pour etre parfaitement modelisee. Des lors que 
l'on restreint le nombre de triangles, on degrade la qualite de son approximation. 
En revanche, sur des regions planes d'un objet, on peut decimer jusqu'a 
l'obtention d'un maillage optimal (ainsi, sur un parallelepipede, l'optimum est 
atteint lorsque chaque face est decrite par deux triangles). 
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On souhaite coder progressivement les maillages, et done retirer 
iterativement les sommets par une operation de fusion d'aretes. Apres 
l'optimisation des positions, un autre processus est done developpe, consistant a 
definir une priorite sur les transformations elementaires, de maniere a effectuer en 
5 . premier les transformations affectant le moins possible le modele. 

Pour cela, on definit la notion de cout pour chaque operation £lementaire, 
et une queue de priorite triee selon ce cout. 

Ainsi,.initialement, toutes les transformations possibles sont listees puis le 
cout de chacune d'elles est calcule (sans effectuer les modifications 
10 correspondantes sur le maillage). On realise la transformation de moindre cout, 
puis on recalcule les couts des transformations elementaires modifiees par la 
transformation precedente. On rappelle que la fusion d'arete supprime deux faces 
soit en consequence, les transformations concernant ces faces sont retirees de la 
queue de priorite, ce qui cree une ou deux nouvelles transformations, pour 
15 lesquelles on calcule egalement le cout en energie, deduit de la technique decrite 
en annexe 2, et qui represente la variation de volume entre le maillage avant et 
apres transformation. 

Le calcul direct du volume entre deux surfaces triangulees etant complexe 
a calculer, on fait appel a une heuristique permettant de simplifier ce probleme. 
20 Soient : 

la transformation fiisionnante de sommet Xj.et Xj de M' ; 
Xi f le sommet de NT resultant de la fusion (dont on a optimise la 
position) ; 

F(X/) les faces de M' avoisinant le sommet X? apres transformation ; 
25 V M Tensemble des sommets de M appartenant aux faces ayant ete 

intersectees lors du calcul de 1' orientation des surfaces pendant 
T optimisation. 

Le cout d'une transformation elementaire T f s'exprime alors sous la forme: 
C(T i (X i ^ i )) = maxd 2 (V M ,F(X i f )) 
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Ce cout correspond done a la distance maximale des sommets du maillage 
original (sous les faces avoisinant le sommet X 4 ) aux faces avoisinant X x apres 
transformation. La figure 16 illustre, en deux dimensions, Texpression du cout 
dans le cas d'une decimation sur une discontinuity. Comme indique 
5 precedemment, en trois dimensions les segments sont remplac& par des triangles. 
La transformation 181 consiste a fusionner les sommets X i et Xj en un sommet 
Xj\ Les sommets 1S2 X a 182 4 correspondent a l'ensemble V M , 

On peut constater que sur un plan maille, le cout d'une transformation 
serai t nul, les zones comportant une information redondante etant decimees en 
10 priorite. 

On presente maintenant des exemples de synoptique pour la mise en oeuvre 
des differents aspects discutes precedemment. 
63.4.7 synoptiques 

• approximation de maillages 
15 La figure 18 resume I'algorithme d'approximation d'un maillage M 

comportant n faces par un maillage NT comportant m faces (avec bien sur m < n). . 

On commence tout d'abord par copier le maillage source M, ou maillage 
original, dans la variable M\ representant le maillage simplifie que Ton cherche a 
determiner (191). Sur ce maillage M\ on effectue ensuite une recherche de toutes 
20 les transformations elementaires realisablcs (192), puis un calcul (193) des couts 
energetiques de ces transformations apres optimisation. 

On selectionne et on realise (194) la transformation de moindre cout parmi 
l'ensemble des transformations possibles. II s'agit de Fetape de decimation. 

Puis, on met a jour (195) les couts energetiques des transformations, dans 
25 le voisinage affect6 par la transformation effectuee a Fetape 194. Enfin, on ajoute 
(196) les nouvelles transformations induites par la transformation effectuee 194, et 
on calcule les couts correspondants. 

Les trois etapes 194, 195 et 196 sont reiterees (197) n - m fois, jusqu'a ce 
qu'au moins un des criteres d'arret soit atteint. 
30 • optimisation en position 
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La figure 20 illustre roptimisation du sommet Xj, apres fusion de l'arete 
■«,Xj). 

On definit au pr£alable les parametres d'initialisation suivants, 
partiellement en fonction de la dynamique de la g^ometrie locale autour du 
sommet X i : 

- Epsilon est un critere d' arret de la convergence, fraction de la longueur 
moyenne des aretes autour du sommet Xj (typiquement 0,001) ; 

- Maxlter est la borne superieure du nombre d' iterations, fixee 
arbitrairement (typiquement 200 iterations) ; 

- Steplnit, StepDec et StepMin sont respectivement la valeur 
d'initialisation, le taux de decroissance et la borne inferieure du pas. 
Steplnit et StepMin sont des fractions de la longueur de l'arete a 
fusionner (typiquement 0,1 et 0,01), StepDec est fixe arbitrairement 
(typiquement 0,95) ; 

- StepSampling est le pas d'echantillonnage sur les triangles, utilise pour 
le calcul discret de 1 'integrate formant le gradient a chaque iteration. II 
est deduit des aires des faces voisines du sommet X. 

- NbSharpEdges(X) est le nombre d'aretes vives autour du sommet X. 
L'algorithme debute par une etape 201 d'initialisation en position, qui 

consiste a effectuer les operations suivantes : 

Si (NbSharpEdge(Xi) = NbSharpEdge^)) 
Alors 

initialiser au milieu de l'arete a fusionner, 

sinon 

initialiser sur le sommet du plus grand nombre d'aretes vives. 
On realise ensuite le calcul 202 du gradient, avec les operations de 
discretisation sur les triangles et d'orientation des surfaces. Puis, on effectue le 
deplacement 203. Le sommet se deplace de la valeur du vecteur gradient multiplie 
par le pas courant. 
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On controle ensuite revolution 204 du pas. Si une oscillation est detectee 
avec le vecteur gradient precedent, le pas est multiplies par StepDec. Si le pas est 
inferieur a StepMin, il est reinitialise a StepMin. 

Enfin, on effectue un test de convergence 205. Si la distance entre la 
position courante et la position precedente est inferieure a Epsilon, ou si le 
nombre d'iterations Maxlter est atteint, l'algorithme est arrets (206). Sinon, on 
r^pete (207) les etapes 202 a 205. 
6.4.3.8 resultats 

On montre maintenant, a partir d'un exemple simple (un parallelepipede), 
le comportement du precede d'approximation de maillages selon 1 'invention. 
L'interet du parallelepipede est de permettre de verifier la qualite de 
Tapproximation sur les discontinues (aretes vives regulieres et coins) et sur les 
plans, 

Sur un tel objet, l'optimum est atteint lorsque Ton obtient chaque face 
maillee par deux triangles et que la forme globale de I'objet est parfaitement 
conservee. 

La figure 19 illustre le cas d'une decimation sur un coin 211 du 
parallelepipede. La solution initiale a ete placee sur le coin 211, puisque le 
nombre d'aretes vives est maximal en ce sommet (trois autour du coin 21 1 et deux 
autour de 1'arete vive reguliere de I'objet (sommet 212)). L'algorithme 
d'optimisation oscille, puis converge autour de cette position initiale 213, puisque 
le pas de gradient a ete initialise a une valeur non nulle, puis a diminue des que 
Foscillation a ete detectee. 

La figure 20 presente une decimation sur une arete vive reguliere. La 
solution initiale a ete placee au milieu de 1'arete 223, puisque le nombre d'aretes 
vives est identique (deux autour de chaque sommet 221 et 222). Plusieurs 
solutions situees sur 1'arete vive reguliere de I'objet satisfont l'invariance de 
volume, et l'algorithme converge apres oscillation autour de 1'arete, sur une 
position dependante de la solution initiale. Le milieu 224 de Tarete constitue done 
une bonne initialisation dans ce cas. 
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La figure 21 illustre une autre situation de decimation sur une arete vive 
reguliere, dans laquelle la solution initiale a ete placee sur l'arete vive 234, 
puisque les nombres d'aretes vives sont differents (deux autour du sommet 233 
situe sur l'arete, et aucune autour du sommet 232 situe sur la face plane 235 de 
l'objet). Plusieurs solutions situ£es sur l'arete vive reguliere de l'objet satisfont 
Tinvariance de volume, et l'algorithme converge apres oscillation autour de 
l'arete 234, sur une position 236 dependant de la solution initiale. 

La figure 22 pr^sente le cas d'une decimation sur un plan 243. La solution 
initiale a ete placee au milieu de V arete 245 formee par les sommets 241 et 242, 
puisque les nombres d'aretes vives sont identiques (aucune autour de chaque 
sommet 24 let 242). Plusieurs sommets situes sur le plan de l'objet satisfont 
Tinvariance de volume. L'algorithme converge, apres oscillation, de part et 
d'autre du plan, sur une position dependant de la solution initiale. Le milieu 244 
de l'arete constitue done, a nouveau une bonne initialisation dans cette situation. 
Le precede decime le maillage initial en utilisant une queue de priorite. 
La figure 23 montre les douze iterations successives de l'algorithme, 
numerotees (b) a ( m), a partir du maillage original numerate (a). 

On constate que les aretes entrainant une invariance de volume sont 
decimees en priorite. L'optimum est atteint pour les sommets demandes sur cet 
objet, puisque chaque face est alors maillee par seulement deux triangles. 

On remarque que les aretes rejoignant deux coins de l'objet ne sont pas 
fusionnees, puisque cette operation aurait entraine une forte degradation de la 
qualite de 1' approximation, les couts et transformations elementaires 
correspondantes etant superieurs aux couts d'une fusion d'aretes sur un plan ou le 
long d'une arete vive reguliere. 
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5 Codage progressif de maillages 

Hugues Hoppe a developpe" une methode fonctionnant par decimation et optimisation des 
points. La decimation est realisee par une transformation elementaire reversible et facile- 
ment codable : la fusion d'arete qui rassemble deux sommets adjacents en on seul sommet 
suppnme deux faces et optimise ensuite la position du sommet resultant de la fusion (figure 

L'auteur a utilise une approche variationnelle permettant la ininimisation du proces- 
sus, la fonctionneUe d'energie associee est definie par un terme d'erreur et un terme de 
jq regularisation : 

E = E error + Erggui 

Ce terme d'erreur caracterise la distance entre les deux maillages M ei M definie lo- 
calement, et garantit la fidelite aux donnees initiales. Le terme de regularisation permet de 
respecter la topologie d'un maillage et d'assurer Tunicate de la solution. Ann d'initialiser 
le probleme, une fois la face a supprimer determinee, trois ininimisations (partant de trois 
positions initiates differentes choisies de maniere a initialiser correctement le processus d'op- 
trmisation) sont realisees (figure 25) afin de retenir la solution correspondant au minimum 
15 d energie de la fonctionneUe. 

Chaque transformation T t supprimant un sommet de la triangulation et modifiant le 
maillage AT,- en M i+1 , engendre une augmentation de l'energie globale qui peut etre calculee. 
Le -cout en energie resultant de cette decimation est le suivant : 

f A&Ti = E Mi+1 -E Mi 

\ LEy. > 0 W 

Au depart toutes les transformations possibles sur le maillage sont listees, puis le cout 
20 energetique pour chacune d'entre elles est calcule. Ces transformations sont rangees dans une 

queue de priorite triee suivant leur cout estime. Au cours de l'algorithme, la transformation 
de cout minimum est effectuee, puis les couts des transformations dans le voisinage affecte 
par cette modification sont mis a jour(figure 26}Les transformations nouvellement creees sont 
ajoutees et leur cout estime est calcule (figure 27) L'algorithme agit done de maniere iterative 
et decime le maillage M d'un sommet et de deux faces a chaque iteration. Cette approche 
de decimation et d 'optimisation est bien adaptee au contexte de codage progressif, puisqu'il 
suffit d enregistrer la sequence des decimations dans le sens inverse (celui des raffinements) 
pour coder un maillage sous une forme progressive. Ainsi un maillage M peut etre represente 
par un maillage a la resolution M 0 et par une sequence de raffinements {r 0) ri,r 2 ,...,r n }. 

L'erreur utilisee pour caractenser la distance entre les deux maillages est la distance 
euclidienne (ou quadratique) entre le point X { a optimiser du maillage M et le maillage M 
qu'U convient de definir. Les points situes a une distance D de Xi dans BP foment une 
sphere de rayon D : 

d2{X u M) = ^d 2 {X i ,Y)=™*\\X i -Y\\ 2 ( 2 ) 
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ou 

N 



\\x\\2 = (E x !) 1/2 >X£® N (3) 

Le maillage M etant une surface triangulee, cela engendre une mesure de la distance entre 
un point X t et un ensemble de triangles. On remarque d'apres la definition de la distance L 2 
dans I 3 que le point de M situe a la distance rninimale de X> n'est pas toujours unique. Le 
choix effectif de ce point depend de I'ordre d'insertion des points et de la nature du test lors 
de la recherche du minimum (< ou <). La metrique L 2 appliquee sur les maiUages devient 
la metrique D 2 puisque Ton raisonne avec une distance point-surface. 
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Optimisation 

La technique d'optimisation choisie pour resoudre notre probleme de minimisation est 
une methode iterative, appelee gradient adaptatif. 

1 Gradient adaptatif 

Notre objectif consiste done a minimiser la fonctionnelle E decrite precedemment. Comme 
dans tous les cas, ia stationnarite de E est une condition necessaire d'optimalite, effective- 
ment, pratiquement toutes les methodes d'optimisation sans contraintes dans E N consistent 
a rechercher un vecteur X stationnaire (VE(X) = 0). Ce probleme est equivalent a la 
resolution du systeme d'equations non lineaires : 



||(*) = 0,Vi = l,..,iV (1) 

On peut chercher a resoudre directement ce systeme, ce qui conduit a la methode de Newton. 
Cependant, cette methode peut ne pas converger, si le point de depart des iterations est trop 
eloigne de X. D 'autre part, elle suppose la fonction deux fois continument differentiable, et 
elle necessite le calcul des derivees secondes en chaque point. 

C'est pourquoi les methodes les plus couramment utilisees proc&dent differemment : il s'agit 
de procedures iteratives ou Ton engendre une suite de vecteurs X° t X 1 , ... convergeant vers 
un optimum local de E. 

Ces methodes de gradient constituent une famille de methodes qui procedent de la facon 
suivante : On part d ? un point X° et on calcule le gradient VE(X°) en X° de la fonction E 
a minimiser. Comme V2?(J£°) indique la direction de la plus grande augmentation de E, on 
se deplace d'une quantite 70 dans la direction opposee au gradient, et on definit le point : 



•_ Y 1-_Y0 T „ ;: Vg ^°) , 9 * 

La procedure est repetee et engendre les points X°, X 1 , X n , ... suivant la relation : 



lk * II V£(Jt* || (3) 

Dans cette famille, il convient de signaler les methodes de gradient a pas determine 
dans lesquelles on choisit a priori les valeurs de deplacements 7*. L 'inconvenient de cette 
procedure est que la convergence peut etre tres lente. Nous avons done choisi d'utiliser la 
methode de gradient a pas adaptatif, qui consiste k dlminuer (resp. augmenter) le pas 7* 
lorsque Terreur augmente (resp. diminue), tout en conservant la contrainte de controle du 
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pas, soit : 7* € [jminilmax]- Finalement, comme entire d 'arret, nous nous sommes fixes que 
la difference entre deux erreurs successives devait etre inferieure a un e donne. 

Au cours de la minimisation, les positions des sommets du maillage vont evoluer de 
maniere a minimiser le volume entre les deux surfaces, cette evolution necessite de calculer 
le gradient a l'iteration k (le terme VE(X k ) de l'equation 



2 Variation elementaire de volume 

Le processus d'optimisation est done vu comme une evolution du maillage M de maniere 
a minimiser F(M,M). L'evolution elementaire du maillage est definie par le champ de 

vecteurs SM. En realisant une parametrisation de la surface par u et r, le vecteur SM 

peut s'exprimer sous la forme SM{u, u). Si on approxime localement la surface par son plan 
tangent la variation elementaire de volume est alors le parallepipede genere par revolution 

de 1'element de surface dudv dans la direction 8M(u,v) (figure 6). 

On peut ainsi exprimer la variation elementaire de volume induite par 1'element d'aire 
dudv : _ 



SV{M,M,SM{u,v)) = 



dM(u,v) dM(u t v) 
du dv 



•Jjfir(u,v) 



• dudv 



La variation est negative lorsque le vecteur SM est oriente vers l'interieur du volume 
d'erreur. la figure 7 presente le cas particulier ou revolution du maillage se fait dans la 
direction de la normale. Prenons maintenant en consideration le fait que : 



dM(u,v) dM(u,v) 
du dv 



• dudv — n(u } v)da(u, v) 



avec n(«, v) la normale unitaire en M(u, v) et da{u y v) = ||^ A ^jf || • dudv Telement d'aire. 
On exprime la variation de distance entre M et M : 

6d{M % M t 6lfc) = J J \ri(u,v)8M(u,v)\dtr(u,v) 

Cette expression n'est pas facilement utilisable a cause de la valeur absolue. Nous intro- 
duisons maintenant la fonction : t](u,v) = t?(M, M,n(u,v)) = 1 si la normale est orientee 
vers le volume interieur et -1 sinon, nous obtenons alors 



d'ou 



|n(u,v)£M(u,t;)| = f}(u i v)n{u 1 v)SM{u,v) 
<Sd(M, M, 6&) = - / f v)n{u, v)S*M{u, v)da{u, v) 
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3 Modele d'elements finis triangulaires 

Afin de calculer le gradient en chaque sommet du maillage a Iteration k de ralgorithme 
d'optimisation, on parametre le maillage par un modele d'elements finis triangulaires P 1 : 



N 



(4) 



Xi est un sommet du maillage defini dans BS 3 et A^u, v) est une fonction de forme adaptee a 
un modele d'elements finis triangulaires. On definit une fonction de forme sur le maillage par 
le biais de coordonnees barycentriques tel que cela est decrit dans le document "la methode 
des elements finis, formulation de base et problemes lineaires" , par O.C. Zienkiewicz et R.L. 
Taylor (afhor technique, 4ieme edition, 1991). Pour cela on numerate les sommets du triangle 
(figure 8), et on definit un systeme de coordormees A : , A 2 et A 3 liees par les relations : 



x — Ai • Xi 4- A2 * X 2 + A3 * 



X3 



y = Ai • y x + A 2 * V2 + A 3 • y z 
1 = A 1 +A 2 + A 3 

Nous allons maintenant determiner les valeurs des fonctions de base : au sommet 1, 
Ai = 1, A 2 = 0 et A 3 = 0. Les lignes de niveau de X x sont des lignes equidistantes paralleles 
au cote 2 - 3 le long duquel \ x = 0. On peut exprimer la coordonnee Aj en un point P par 
un rapport d'aires : 

_ aire{P2Z) 
~ aire(l23) 



Ai = 



(5) 



On obtient 



Ai = 
A 2 = 
A 3 = 



ai + bi • x -f c\ • y 





2A 


a 2 4- 6 2 


•z + c 2 -y 




2A 


a.3 + &3 


•x + cz-y 



Avec : 



2A 



A = azVe(123) = - • det 



Xi 
X2 



V2 



1 S3 y 3 
<*i = 22 • Vz - x 3 • y 2 
h = y 2 - 2/3 
Ci = X3 — x 2 

On calcule les valeurs a 2 , 62, c 2 , a 3 , 63, et c 3 par permutation circulaire des indices 1, 2 et 3. 
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4 Evolution de surface 

L'evolution du maillage peut ainsi s'ecrire : 

i=i J Ju > v 

on en deduit la derivee partielle de la distance pour un sommet X { : 
dV _ f f 

L'algorithme de gradient pour un sommet X { s'ecrit sous la forme : 

ou 7i represente un pas adaptatif avec 7 i 6 [7mm, To]* Au depart le pas est initialise^ a 7o 
P ^ + f t m ^ iipM par k ( avec k < 1 ) en Presence d'une oscillation, c'est a dire lorsque 
< "axT' Wn: ><0Le comportement de l'algorithme d'optimisation peut-etre illustre en 2d 
(figures 10 et 11), revolution du sommet X { du modele M est liee a la somme des vecteurs 
d'orientation du modele par rapport a la courbe originate M, ponderes par la fonction de 
forme qui vaut 1 en X { et decroit en 0 le long des segments voisins de X it En 3d le principe 
est le meme avec une fonction de forme definie sur les triangles voisins du sommet X>. 
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REVENDICATIONS 

1. Procede de simplification d'un maillage source M forme d'une pluralite de 
surfaces d^finies par des sommets, des faces et des orientations sur ces dernieres, 
ledit procede mettant en oeuvre une etape de decimation par fusion d'arete, 
consistant k associer & une arete a decimer, definie par deux sommets (1 1 , 12), un 
sommet unique ( 1 3) , de fafon a obtenir un maillage simplify M\ 

caract&ise en ce qu ! il comprend une etape de pseudo-optimisation apres ladite 
etape. de decimation par fusion d'une aretes, positionnant le sommet resultant de 
ladite fusion de maniere a reduire la deviation geom6trique entre ledit maillage 
source M et ledit maillage simplifie M 1 . 

2. Procede de simplification d'un maillage source selon la revendication 1, 
caracterise en ce que ladite etape de pseudo-optimisation consiste a denombrer les 
aretes vives autour des deux sommets formant Tarete a fusionner, et a distinguer 
les deux cas suivants : 

si les nombres d'aretes vives sont les memes autour des deux 
sommets, on place le sommet resultant de la fusion au milieu du 
segment reliant lesdits sommets (42, 44) ; 

si les nombres d'aretes vives sont differents, on place le sommet 
resultant de la fusion sur le sommet presentant le plus grand 
nombre d'aretes vives (43, 45). 

3. Procede de simplification d ! un maillage source selon Tune quelconque des 
revindications 1 et 2, caracterise en ce qu'il comprend une etape de selection 
d'une fusion d'arete a effectuer, parmi toutes les fusions d'arete possibles, tenant 
compte : 

d ! au moins une information representative de la courbure definie 
localement autour de l'arete consideree ; 

d'au moins une information representative de la dynamique 
geometrique definie localement. 

4. Procede de simplification d'un maillage source selon la revendication 3, 
caracterise en ce que ladite etape de selection met en oeuvre une queue de priorite 
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des aretes a fusionner, en fonction d'un critere prioritaire, ladite information 
representative de la courbure, puis d'un critere secondaire, ladite information 
representative de la dynainique geometrique. 

5. Procede de simplification d'un maillage source selon la revendication 4, 
5 caracterise en ce que ladite etape de selection gere un seuil de courbure, seules les 
aretes ayant une courbure inferieure audit seuil etant considerees pour 
l'application dudit critere secondaire, 

ledit seuil etant augmente lorsqu'aucune arete ne presente plus une courbure 
inferieure a ce dernier. 

10 6. Procede de simplification d'un maillage source selon Tune quelconque des 
revendications 1 a 5, caracterise en ce que ladite information representative de la 
dynamique geometrique appartient au groupe comprenant : 
longueur de l'arete consideree ; 

une moyenne des surfaces des faces avoisinant ladite arete 
15 considerfe ; 

une moyenne des longueurs des aretes adjacentes aux sommets 
formant ladite arete consideree ; 

une combinaison de longueurs d'aretes et/ou de surfaces de faces. 

7. Procede de simplification d'un maillage source selon Tune quelconque des 
20 revendications 1 a 6, caracterise en ce que la decimation est interrompue en 

fonction d'un des criteres appartenant au groupe comprenant : 
un taux de compression atteint ; 

une complexity geometrique atteinte, exprimee par un nombre de 
sommets ou de faces ; 
25 - un seuil de courbure atteint. 

8. Procede de simplification d'un maillage source selon l'une quelconque des 
revendications 1 a 7, caracterise en ce qu'il constitue une etape d' initialisation 
d'un procede d' optimisation geometrique d'un maillage. 
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9. Procede d'optimisation geometrique d'un maillage source, caracterise en 
ce qu'il comprend une etape d' initialisation mettant en oeuvre le procede de 
simplification de Tune quelconque des revendications 1 a 7. 

10. Proc&te de codage d'un maillage source (M) selon la revendication 9, 
representatif d'un objet en trois dimensions, delivrant un maillage simplifie (M') 
correspondant audit maillage source (M), lesdits maillages etant definis par un 
ensemble de sommets, d'aretes et/ou de faces, 

caracterise en qu'il met en oeuvre une etape de simplification selon Tune 
quelconque des revendications 1 a 7, puis une etape de minimisation d'un volume 
compris entre ledit maillage source (M) et ledit maillage simplifie (M'). 

11. Procede de codage d'un maillage source selon la revendication 10, 
caracterise en ce que, chacun desdits maillages etant defmi par la position de 
chacun de ses sommets, ladite etape de minimisation assure la determination des 
positions des sommets -(X,, X 2 ,... X n ) dudit maillage simplifie (M') minimisant le 
volume V(M, M') entre ledit maillage source et ledit maillage simplifie. 

12. Procede de codage d'un maillage source selon Tune quelconque des 
revendications 10 et 11, caracterise en ce que ladite etape de minimisation met en 
oeuvre un processus iteratif optimisant progressivement les positions des sommets 
dudit maillage simplifie (M'). 

13. Procede de codage d'un maillage source selon la revendication 12, 
caracterise en ce que ledit processus iteratif est interrompulorsqu'au moins un des 
criteres d'arret suivants est atteint : 

un nombre maximum d' iterations ; 

une difference entre deux vecteurs deplacement successifs des 
positions des sommets inferieure a un seuil predetermine (Epsilon). 

14. Procede de codage d'un maillage source selon Tune quelconque des 
revendications 1 1 a 13, caracterise en ce que ladite etape de minimisation met en 
oeuvre une methode de gradient adaptatif . 
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15. Precede de codage d'un maillage source selon la revendication 14, 
caracterise en ce que ladite methode de gradient adaptatif repose sur les operations 
suivantes: 

selection d'un vecteur Xpde R 3n (n s 1) dudit maillage simplifie, et 
5 calcul du gradient VE (Xp) en Xp de la fonction a raininiiser : 

E=d(M,M'(X„...^ I1 )); 
determination de la position Xp" de Xp dudit maillage, selon la 
relation definie a l'iteration k+ 1 par : 

Xp =Xp " YkX Fix7^| 

10 k variant de 0 a n-l etY k etant le pas de ladite relation. 

16. Procede de codage d'un maillage source selon Tune quelconque des 
revendications 10 a 15, caracterise en ce que, a chaque iteration, on determine une 
variation elementaire dudit volume correspondant a un champ de vecteurs ^ , 

et en ce que, la surface etant parametree par u et v, de fa$on qu'un vecteur 6 ^ 
15 s'exprime sous la forme 5M * (u, v), ladite variation elementaire est assimilee au 
parallelepipede genere par revolution de 1'element de surface dudv dans la 
direction ^ (u, v). 

17. Procede de codage d'un maillage source selon Tune quelconque des 
revendications 10 a 16, caracterise en ce que ledit maillage simplifie est param&re 

20 a Taide d'un modele d'elements finis. 

18. Procede de codage d'un maillage source selon la revendication 17, 
caracterise en ce que lesdits elements finis sont obtenus a 1'aide d'un interpolateur 
affine. 

19. Procede de codage d'un maillage source selon l'une quelconque des 
25 revendications 10 a 18, caracterise en ce qu'il met en ceuvre un codage progressif 

dudit maillage simplifie, par decimation d'un sommet et optimisation locale ou 
par decimation d'un ensemble sommets et optimisation globale. 

20. Procede de codage d'un maillage source selon Tune quelconque des 
revendications 10 a 19, caracterise en ce qu'il comprend une etape de limitation 
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de la degradation due a une transformation elementaire, mettant en oeuvre une 
queue de priori te sur les transformations elementaires. 

21. Precede de codage d'un maillage source selon la revendication 20, 
caracterise en ce que ladite etape de limitation de la degradation due a une 
transformation elementaire, e'est-a-dire une fusion d'arete, definie par deux 
sommets, comprend les e tapes de : 

calcul d'un cout pour chaque transformation elementaire possible ; 
realisation de la transformation elementaire de moindre cout ; 
recalcul des cofits des transformations elementaires modifiees par la 
transformation elementaire precedente ; 

ajout des nouvelles transformations elementaires creees et calcul des 
couts correspondants. 

22. Procede de codage d'un maillage source selon la revendication 21, 
caracterise en ce que le cout d'une transformation elementaire (T { ) s'exprime par : 

CO\(X . ,X j)) = max d 2 (V M ,F(X f)) 

avec : T { transformation fusionnant deux sommets X, et X } dudit 
maillage simplifie M' ; 
X/ le sommet dudit maillage simplifie M' resultant de ladite 
transformation ; 

F(Xi f )les faces dudit maillage simplifie M' avoisinant le sommet X/ 

apres ladite transformation ; 
V M ensemble des sommets dudit maillage source M appartenant 

aux faces ay ant ete intersectees lors du calcul de 1' orientation 

des surfaces pendant ladite minimisation. 

23. Application du procede de codage d'un maillage source selon l'une 
quelconque des revendications 1 a 22 a au moins un des domaines suivant : 

la realite virtuelle ; 

la simulation scientifique ; 

la modelisation. 
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